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In this master thesis, a new approximation scheme to non-relativistic potential scattering 
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tions of the T-matrix, which permit the application of the Feynman-Jensen variational 
method. A simple Ansatz for the trial action is made, and, in both cases, the varia- 
J> tional procedure singles out a particular one-particle classical equation of motion, given 

in integral form. While the first is real, in the second representation this trajectory is 
complex and evolves according to an effective, time dependent potential. Using a cumu- 
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The high energy behavior of the approximation is investigated, and is shown to con- 
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higher terms. Our results are then numerically tested in two particular situations where 
others approximations turned out to be unsatisfactory. Substantial improvements are 
found. 
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1. Introduction 

The importance of variational methods in physics cannot be overestimated. The variational 
principle itself, which states that laws of nature are such that an appropriate function, or 
functional, is extremalized, can be safely traced back at least to ancient Greece. A that 
time. Hero of Alexandria formulated that the path of a light ray in a medium is the one that 
minimizes its length, and derived, from these considerations only, the reflection law in optics. 
In the 18th century, after works of Euler, Lagrange and others, the principle of least action 
was born. The principle states that the physical laws can be deduced from the condition of 
stationarity of a single quantity, called the action S, 



This approach proved so fruitful that nowadays, in most of modern physics, the starting point 
of almost any theory is an appropriate action, and the dynamical laws of the theory are 
determined through the requirement that the action is extremalized. 

Unfortunately, very few physical problems can be solved exactly, but variational methods 
turn out to be very useful as approximation procedures, too. Certainly, there exists as many 
ways to approximate solutions to a given problem as problems themselves, but variational 
approximations are conceptually a class of their own, and can furnish fairly reliable results. 
Let us imagine that we are interested in a function S[x], at a precise point Xph that represents 
schematically the physical system of interest, but, for some reason, our knowledge of the system 
is incomplete, and we do not know what x^^ is. What we know is only that the quantity S[x] 
is stationary at this particular point Xp^^. The very idea of this approximation is to restrict 
the range of values that can take x, namely, we require x to lie on some particular subspace of 
the space in which x took its value initially. Let us parametrize this subspace with the help of 
some parameters a,b,... so that x = x{a,b, ...). If the subspace in question is sufficiently broad, 
the point Xt at which S is stationary with respect to the parameters a, b,... may furnish a good 
approximation of Xph, and, what is more important, an even better approximation of S'[xph], 
since this function is stationary a this point. We have thus reduced the task of finding Xph 
to a simpler optimization problem. Nevertheless, success is not guaranteed, since it depends 
crucially on the chosen subspace, which may or may not be sufficiently close to the physical 
configuration. 

Perhaps the most well-known variational approximation in quantum mechanics is the Ritz 
method to obtain approximate energy levels of bound quantum systems. An account can be 
found for example in [1] . To emulate the unknown wave function of a quantum system, a trial 
wave function depending on some parameters is introduced. The true ground state energy 
£"0 of the system can be shown to be always lower than the expectation value of the system's 
Hamiltonian H, 



The optimization procedure then furnishes an approximative eigenstate, and a value for the 
corresponding energy that is even closer to the exact one. The fact that the true ground 
state energy always lies below the approximative value obtained, makes this method especially 
useful. More than a variational principle, is is a minimum principle. 



6S = 0. 
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In this master thesis, we apply such a variational method in the framework of poten- 
tial scattering, in non-relativistic quantum mechanics. This is an old, time-honored subject 
whose treatment can be found in many textbooks [2], [3]. Potential scattering theory con- 
siders the elastic scattering of a quantum mechanical particle, that interacts through a local 
potential, which itself falls off sufficiently rapidly, so as to permit the existence of free initial 
(long before the scattering) and final (long after the scattering) states. The objects of interest 
are the scattering matrix (S-matrix), that expresses the probabilities of a transition between 
given initial and final states, and (or) the closely related T-matrix. Potential scattering theory 
has the advantage of being numerically solvable through the partial-wave expansion. We use 
therefore this framework in this thesis as a test case for more complicated processes, such as 
many-body scattering. 

Approximation methods in potential scattering theory are numerous. To mention some, a 
truncated partial-wave series is very useful for low energies, while the eikonal expansion[4] or 
the Born expansion cover the high-energy range. It has to be said that variational methods 
in potential scattering exist already, too. One could cite for example a method introduced by 
Kohn [H] in 1948, that is now abundantly used in electron or positron scattering from atoms 
or molecules. Also, a stationary expression for the T-matrix was given by Schwinger p], which 
can be used to calculate the partial-wave phase shifts. So, what brings this thesis that is new? 

The starting point are two exact path-integral representations of the T-matrix, that were 
very recently published [7j. The path- integral formulation of quantum mechanics was first 
worked out in the 1940's by Richard Feynman [S], based upon a paper by Paul Dirac on the 
role of the Lagrangian in quantum mechanics. This formulation brought new insights into the 
structure of quantum physics, and proved essential to the subsequent development of theoret- 
ical physics. In this formulation, the probability of an event is explicitly written as a sum over 
all paths, that is over each and every possibility that leads to that event. Each path carries a 
probability amplitude, a complex number, that has the same magnitude for each path, but has 
a different phase, given by i times the ratio of the action S to the reduced Planck's constant 
h. This formulation naturally contains classical mechanics, in the sense that if the action is 
much larger than h, the paths that contribute most to the realization of the event are those 
for which the action is stationary, i.e. the classical paths. Besides, many applications were 
subsequently found in statistical physics or quantum field theories, where it is now a major tool 
|9][T0]. Path integrals are convenient in the sense that they deal with classical quantities and 
usual numbers instead of quantum mechanical operators. A variational principle for actions 
in a path integral formalism is provided by the Jensen inequality for convex functions. With 
the help of a trial action 5*^, one shows that 

where the average is taken with respect to the positive, convex, weight function e~^*^^, and 
AS stands for the difference S — St. This leads, as in the case of the Ritz method, to a 
minimum principle, that is used for instance in statistical mechanics, and was first introduced 
by Feynman. The main restriction to the method is the fact the trial action must be simple 
enough to allow analytical computations to be pursued. Corrections to this approximation 
can also be evaluated. 
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The two path-integral representations mentioned above for the T-matrix are free of any 
seemingly diverging phases in the limit of asymptotic times, and instead avoid the need of 
infinite limits that occurred in previous descriptions, which rendered identities rather formal 
and not well suited for practical calculations. While it was shown in that these represen- 
tations lead immediately to a new high-energy expansion scheme, they have the virtue that a 
variational procedure is also possible. This particular variational approximation is very sim- 
ilar to the one we just described above, and has the valuable advantage over the previous 
variational approximations we quoted, that it is not tied to the partial wave formalism, nor 
is restricted to high energies only. Also, it is the action that is forced upon the optimization 
procedure, and not a wave function, or a T-matrix element; it forms thus a completely new 
way to address the scattering problem, which permits to gain new insights. The difference 
with the Feynman- Jensen procedure is, as we will explain in more a detailed fashion later, 
that we will deal with complex quantities, i.e scattering amplitudes. As a consequence, our 
variational principle will not be a minimum principle anymore. This is a property also shared 
by the other quantum mechanical variational principles of Kohn and Schwinger. 

The structure of the thesis is the following. In section 2, the stage is set up: we derive 
the two exact path-integral representations of the T-matrix. The derivation follows closely [7] . 
In section 3, a stationary expression for the T-matrix is written down, the Feynman- Jensen 
variational principle and our variational Ansatz are presented in some more details. Correc- 
tions are also discussed. In section 4 and 5, the variational approximations to the T-matrix 
corresponding to the two representations are performed and discussed. The first correction 
to the approximation is also given. Numerical tests for two potentials are then presented in 
section 6, followed by our conclusions and outlook. Some technical computations are left for 
the appendices. 

In this work, we use natural units in which h = 1. 
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2. Path Integrals for the T- Matrix 



In the framework of nonrelativistic potential scattering, we consider a potential V^(r), that 
vanishes at infinity, initial and final free states (pi and 0/, to which is associated an initial 
and final momentum kj and kj. The S-matrix is the matrix element, taken between these 
free states, of the quantum mechanical propagator in the interaction picture, evaluated at 
asymptotic times. 



5. 



lim 

T-5>oo 



Ui{T,-T) 



(2.1) 



The free states are normalized through {(j)f\(j)i) = (27r)^5(kj — kj). The definition of the 
propagator is 

(2.2) 



(2.3) 



From the S-matrix one usually substracts the identity, that does not represent actually the 
scattering process, and factors out an energy conserving Dirac delta function. What is left is 
called the T-matrix: 



U {tb, to) = exp ( -iH{tb - ta) 
while in the interaction picture it reads 

Uiitb, ta) = exp (iHotb ) exp ( -iH{tb - ta) ) exp ( iHota 



Si 



i2nf6iK - k/) - 27r tS{E, - Ef)Ti^f. 



(2.4) 



We can now start the derivation of the path integral formulation of the T-matrix. We will follow 
closely |7], and present the main points of the derivation. We refer to it for a more complete 
discussion. The starting point is the Feynman-Kac formula, that presents the matrix elements 



U {tb, Xb, ta, Xa) := (x;, exp ( -iH{tb - ta 



(2.5) 



of the propagator between eigenstates of the x operator as a functional integration over paths. 
It formally reads 



U (tb, Xfo, ta, Xa) 



3N 

( m \ — 




(2.6) 
(2.7) 



In this expression xq is meant to be x^, and xat is Xf,. This identity relies on the subdivision of 
the time intervall {tb — ta) into N small steps of duration e = {tb — ta)/N , and the insertion of 
a complete set of x eigenstates between each such step. One thus integrates over all possible 
intermediate values between x^ and x;,. We will however use a different version of this formula, 
i.e. one that uses integration over velocities, and not paths. 
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Velocity Path Integrals This is done by inserting a "big one" in the Xj integrations, in 
the form of 

N 



i=Yijdv,s(v,-'^'^^P^y {2.1 



The (5-functions so introduced enforce sucessively 

3 

Xj=Xo + e^Vj, (2.9) 

i=l 

and permit to effectuate each Xj integration. However, from the N 5-functions, one is still left, 
that enforces 

Xb = x„ + e2^Vj, or x^ = x„ + / div(i), (2.10) 

i=l 

in continous notation. One has therefore 

U{tb,Xb,ta,Xa) = Jim^(^^--^j' / Civi---Civjv5 I Xfe-Xa-e J^Vj 

\ i=i 

■ exp(«f:||vj-y(x„ + .^v.)j). (2.11) 

Prom now on we will exclusively employ the continous version of formulae of this type, that 
is, in this case, 

C/(ib, Xb, ta, X„) = J Vv S (^b-^a- dtv{t)^ 

■ e^pl^tj^t {'^v'{t)-V{^{t))]y (2.12) 
Here the T>v term is defined in the discrete version as 

^ 3JV N 



fc=i 



It has the property that the pure gaussian functional integral is normalized to unity. The 
argument x{t) of the potential is 

x(i) = x„ + r dt v{t) = + 1 r ds sgn(t - s)v{s). (2.14) 

This identity can be inserted into the S-matrix. Indeed, after recognizing that 

Si^f = hm e<^'+^^)^ Uf U{T, -T) 0,) , (2.15) 
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one inserts two complete sets of eigenstates and x;, of the x operator, and recovers the 
matrix element ?7(tfe, x;,, x^). Together with (2.12), it results in 

Si^f = ^lim e<^'+^^)^ j d^a j d^b exp {-i [kf ■ X5 + • xj) (2.16) 
• J V^5(^,-^a- dt^{t)^e^v(^ f^dt |%2^t)-\/(x(t))}y 

After the change of variables r = (x + y)/2 and s = x — y, the relative coordinate s is fixed by 
the 5-function. We furthermore define the mean momentum and momentum transfer through 

K = ^(k,, + k;), q = k;-ki, (2.17) 
and obtain the representation 

S,^j = Ymv e<^"+^/)^y" dre-'''-' J exp j^^ dt yv2(t)j 

■exp(^-ij K- v(t) + V^(r + x,(t))^ , (2.18) 



where we used 



1 



Mt) = 2 y sgn(t - s)v(s). (2.19) 

With a simple shift one gets rid of the linear term K ■ v(t). Indeed, with \(t) — )■ v(t) + K/m, 
recognizing that Ei + Ef — K^/m = q^/4m, and that 

r+T 

J dssgn{t- s) = 2t, te[-T,T], (2.20) 

one has finally 

■exp(-z/ V (r + -t + x,(t) ) ) . (2.21) 



-T 



m 

For a vanishing potential one calculates from this expression 

Si^f = (271)' 6iki-kf), (2.22) 
so the part of the S-matrix that describes the scattering is 



■|exp^-iy V (r+^t + x„(t)]^ -iV (2,23) 
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Diverging Phase Is it possible to get rid of the ill-defined phase Tq^/ (4m) in front of the 
scattering matrix? To this aim, the factor is transformed into a 3-dimensional Laplacian 
A acting on the variable r, in the exponent of exp (— iq ■ r). Since we are assuming that the 
potential vanishes at infinity, we can then integrate by parts without any boundary terms and 
let it act on the potential term. One can further reduce it to a shift operator with the help of 




= Jvwexp(-ij dtjw^{t)±j ^/(t)w(t) ■ vV (2.24) 



_rp Z J_rp 

In this expression, f{t) can be any function, as long as it fulfills the condition 



dtf{t) = 2T, 



(2.25) 



so as to recover T times the Laplacian, and the functional integration is normalised as before, 
that is, the pure gaussian integral is the identity. The additional minus sign in the quadratic 
term is mandatory in order to get a real shift operator acting on the potential, while the linear 
term can have any sign. For reasons of convenience we choose the minus sign, and f{t) to be 
sgn(— t). In that case the S-matrix becomes 

{S - 1),^/ := ^lim j dr e"*'!" j V^^ J Vw exp (^t ^ dt ^ [v2(t) - w\t)]^ 

exp(-t r y fr + -t + x,(t)-x^(0)^Vll- (2-26) 



T 



m 



These new degrees of freedom, that appear with a kinetic term of the wrong sign, thus remove 
all infinities. The 3-dimensional vector w will be called antivelocity on some occasions. 

Extraction of the T-Matrix In order to reach for the T-matrix, one needs to extract the 



energy conserving (5-function from (2.26). This is done first by observing that the action in 



the path integral is (at least formally) invariant under the transformation 

t = t + T, r = r-— r, and v(t) = v(f). (2.27) 
Indeed, such a transformation only affects the integration bounds, 

dtV {r + —t + ^^{t) -^^{0) 



-T 



m 



T-T / K- 1 /■+^"^ - \ 

dtV[r + —t+- dsv(s)sgn(t-s)-x^(0) . (2.28) 

-T-T \ m 2 J_T_^ J 

We expect that in the infinite T limit, this change of the bounds does not modify the action. 
This also means that it does not depend on the component of r that is parallel to K. The 
integration over it thus produces a singularity, which is the energy conserving 5-function we 
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are looking for. To extract it, we use the Faddev-Popov trick, i.e. we insert first the "gauge 
fixing" term 

!• / r n \ 

(2.29) 





K 


\ = — f dTdlK- 


r H r 


rn J \ 


m 



in the r integration of (2.26). The vector K = K/if is the unit vector in the K direction. We 
then perform the transformation (2.27), which, as we argued, does not change the action, but 
only the exp(— iq ■ r) in the r integration. The only dependence on r is now to be found in 

K 



dr / dr exp — iq ■ r + i q - — r 5 ( K ■ r 



m 



{2n)6 



q K 



m 



d% exp (— zq • b) 



(2.30) 



As desired, the argument of the (5-function is the difference in energy before and after the 
scattering. The 2-dimensional vector b (the impact parameter) is the part of r that is perpen- 
dicular to K. As a consequence of the conservation of energy, the relations 



K = 2k cos 
hold, where 6 is the scattering angle. 



q = ksin 



(2.31) 



Eikonal Representation of the T-Matrix We have thus found our first representation 
of the T-matrix. It now reads 



n(3-3) 



K 



m 



j d^e-'''-^ J VwVw exp (^i^ j dt [w^{t)-w{tf] 



exp 



dtv{u{t)) 



(2.32) 



where the trajectory is 



U{t) :=b + -t + x„(t)-x^(0). 

m 



(2.33) 



The upperscript (3-3) indicates the presence of the two functional integrations over the 3- 
dimensional variables v and w. We will see that it is possible to write a similar formula with 
an 1-dimensional antivelocity w only. 

The structure of the T-matrix is transparent: the reference trajectory along which phases 
are collected is the straight line, eikonal trajectory b -|- t K/m, while the functional integrals 
describe the quantum fluctuations around this path. We will use many times the notation 



XK(b,v,w) :=- / dtV{U{t)) 



(2.34) 
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Ray Representation of the T-Matrix The fuU 3-dimensional path integral over w was 
introduced with the hope to remove all infinities. We will see in the subsequent chapters that 
this works in a very efficient way. However, it is also possible to use a one-dimensional, longi- 
tudinal, functional integration only, which has the effect of changing the reference trajectory. 
Let us consider in (2.32) the following shifts of variables: 

v(t) = ^sgn(t)+v(t), w(t) = ^sgn(t)+w(t), (2.35) 
and b = b - Xi)x(O) + x^x(O). (2.36) 
A small calculation shows that, on one hand, 



m 
~2 



/TTl f 
dt [v2(t) - w\t)] = ^ I dt [v2(t) - w\t)] + q ■ [x^^(O) - x,x(0)] • (2.37) 



This is to say that the shift (2.36 ) of the impact parameter exactly compensates for the 
difference and we recover the same free action. On the other hand, using the relation 

r+T 

J ds sgn{t - s)sgn{u - s) = 2T - 2\t - u\, t,ue[-T,T], (2.38) 
the argument of the potential, that we will call now Crayit), becomes 

Uyit) = h+-t + ^\t\+ x,^(t) - x,^(0) - x^//(0). (2.39) 

TTh ^TTh 

The dependence on the components of w that are perpendicular to the mean momentum K 
is thus to be found only in the free part of the action. Given our normalisation of Pw, we can 
simply forget about them. The new reference path is now 

b + -t + := b + "^I^t = b + k,e(-t) + k.e(t). (2.40) 

m 2m m 

A particle on this trajectory follows a straight line, according to its initial free motion, suffers, 
at b and at time t = 0, an elastic bump and moves again freely, in the direction given by its 
final momentum kj. This path may intuitively represent better the picture of the scattering 
process than the straight-line eikonal path we have just seen. Besides, it has the advantage that 
the particle has the right scattering energy fc^/ (2m) instead of the slightly misplaced K'^ / (2m) 
in the eikonal version. We call this representation of the T-matrix the ray representation. It 
reads 

tI^P = d'^be-'^ '' J Vw j Vw exp (^"^ J dt y{t)-w{ty] 

■{exp[ixray (b,v,w)] - 1}, (2.41) 



with the phase 



Xray{h, \,w) = - j dt V {^ray 



it)) . (2.42) 
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Impact Parameter Representation It has to be emphasized that both of these represen- 
tation for the T-matrix are not impact parameter (Fourier) representations in the strict sense, 
since the integrand has a dependence on both b and q, the latter through 



At this point, it is appropriate to mention the eikonal expansion of the T-matrix by Wallace 
[4j, which is a systematic high energy expansion of the scattering process, and possesses some 
properties analogous to our approach. For instance, up to second order in inverse powers of 
the incoming momenta k, it reads 



where the subscript accompanying each term indicates the inverse power of k attached to it 
(while the velocity v = k/m is treated as a constant). Similarly as in our representations, the 
eikonal approximation consists of integrating so-called scattering phases over a two-dimensional 
impact parameter vector b. However, unlike ours, the eikonal representation is indeed an 
impact parameter representation, and the factor in front of the T-matrix is k/m and not 
K/m, a difference that plays a role for large scattering angles. Since we will compare our 
own results to the eikonal phases later, let us mention them now, for a spherically symmetric 
potential, (general expressions can be found in or later in this work, in another form) 




(2.43) 




(fh e-'^ '' {exp (z [xo(b) + nih) + r^ih)] - u^ih)) - 1} 



(2.44) 




(2.45) 



(2.46) 



1 1 

2 ^2 ^3 




/ 

^0 



oo 




24A;2 



(2.47) 





where Aj, is the two-dimensional Laplacian acting on b. 
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3. The Feynman- Jensen Variational Principle 



The two formulae ( |2.32[ ) and ( |2.41 ) contain a lot of physics. On one hand, it is possible to 



show (see [7]) that, in both representations, a developpment in series of the exponential term 
that contains the potential leads to the Born expansion, term per term. On the other hand, 
in [7] an high energy approximation of the T-matrix is also done, by expanding the phases 
Xk, or Xray, simultaneously in powers of v and w. This leads to the eikonal expansion we just 
discussed, and eikonal-like approximations. However, and this is the topic of the thesis, these 
identities can also serve in a third kind as an approximation procedure, namely, one can make 
use of the Feynman- Jensen variational principle. This principle is based on the fact that, for 
a given action 5* that depends on x, the quantity 



Vxe-^ (3.1) 



can be exactly rewritten as 



Vx e'^ = ? / Vx e"-^* 



/ Vx e-^* 

^g-(S-5t)^ J jy^^-St^ (3^2) 



for any other action St, which is, in practice, chosen so that one can solve the path integral on 
the right side of the last term. One can use then the Jensen inequality for convex functions 
on the term on the left side, 

(exp [-{S - St)]) > exp [- {S - St)] , (3.3) 



so that one gains an approximation of the quantity (3.1), together with a lower bound, since 
the terms that enter can be readily calculated, if the action Ansatz is suitable enough. Such 
an approach is used, for example, in statistical mechanics, where path integrals in euclidean 
time show up. 

However, our situation does not correspond exactly to this formulation. The quantities that 



we have in (|2.32|) and (|2.41|) are of the form 

Vx e'^, (3.4) 



with (for the eikonal representation), 

Vx = PvPw, (3.5) 

and 

S='^Jdt (v2(t)-w2(t))+XK(b,v,w). (3.6) 



In that case it still holds 



Vxe'" = ^ / Vxe 

J Vx e'^^ 
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However, the weight function over which is taken the expectation value {|e**^'^~'^*^) is now 
a complex function , so one cannot use the Jensen inequality, and looses the lower bound 
information. Nevertheless, although it is not a minimum, the corresponding quantity still is 
stationary with respect to St: Let us define the functional 



F[St] 



Obviously, 



exp(z (S-St)) j Vxe'^' 
exp {i {S - St)) mo 

F[S]= [ Vx e'^ 



{3.i 



(3.9) 



is the quantity we are looking for. To prove stationarity, let us consider a differential change 
in St given by SSt. Then one has that 



6F\ 



St=S 



5mo + imo 5 {S - St)\ 



St=S ■ 



It holds further, 



5mn = 6 I Vxe'^' =i I Vx e'^'6St 



while 



S{S-St)\ 



6- 



. . - . I Vxe'^HS-St] 
st=s I y 



S=St 



Vx e'^'SSt 



so that one effectively has in (3.10) 



(3.10) 
(3.11) 
(3.12) 
(3.13) 



We have thus proven the stationarity of (3.8) at S. Therefore, although one does not find a 



minimum or a maximum, as we said, one can hopefully still have a good approximation by a 
suitable choice of action St. 



Our Ansatz The trial action St has to be chosen so that the quantity mo and the average 
{S — St) can be exactly solved. This restricts the possibilities in the class of quadratic functions 
of their arguments. As we have mentioned at the beginning of this section, an high-energy 
expansion follows from expanding the scattering phases xk and Xray in powers of v and w. 
The effect of such an approach to first order is to render the interacting part of the action 
linear in v and w. We will investigate in this work a simple Ansatz for the action, which is to 
add a linear term in v and w to the free term of 5* that we have. Our approach is therefore 
similar to the high energy expansion, except for the fact that the optimization procedure will 
select for us not the linear term corresponding to this high energy expansion, but the best 
possible one. We expect thus our approach to have a greater validity than for high energies 
only. The trial action therefore reads 



5i(v,w) = ^ j dt (v2(t) -w2(t)) + j dtB{t)-Y{t) + j dtC{t)-w{t) 



(3.14) 
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in the eikonal representation, while in the ray representation, 



St{v,w) 



m 



j dt {^r\t) - w\t)) + j dt B{t) ■ v(t) + J dt C(t) ■ w{t). 



The B's and C's are functions that will be determined by the stationarity of F: 



6F 6F 

0, TT^=0, 



6B{t) 



6C{t) 



respectively 



5F 5F 

0, = 0, 



6B{t) ' 6C{t) 
and the T-matrices will then be approximated by 

.K 



(3-3) 



j d%e-'^ '' (F[5i(v,w)]-1), 

/ d%e-'^'° (F[5i(v,u;)]-1). 
m J 

The corresponding expectation values that enter F and need to be computed are 



is: 



{{S-S,)) = (xk)-((B,v))-((C,w)) 
= (x™,)-((B,v))-((C,t.)), 

where we introduced the short-hand notation 



that we will use frequently from now on. 



dtf{t)-g{t), 



(3.15) 



(3.16) 
(3.17) 

(3.18) 
(3.19) 



(3.20) 
(3.21) 



(3.22) 



Corrections Corrections to this approximation can be gained by recognizing that the quan- 
tity 

exp [i {S - St)] 
is the first term in an expansion in cumulants of 



entering equation (3.7). The cumulants of the distribution are defined [TT] through 



^git(5-50^ = exp 



On the other hand, the expansion in moments 



,fc=l 



(3.23) 



fc=0 



(3.24) 
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permits, through comparison of powers of to express the cumulants in terms of the moments. 
One has for the first two, 

Ai = {(S-St)) (3.25) 
A2 = {{S - Stf) - {{S - St)f . (3.26) 

In this work we will in due time consider the first correction A2, with the corresponding 
functional 

F[St] = mo exp {t {S - S^)) exp (-^ [{{S - S,)') - {{S - St)f]\ . (3.27) 
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4. Variational Approximation in the Eikonal Represen- 
tation 

In this section, we apply the Feynman- Jensen variational principle to the eikonal representation 



(2.32) of the T-matrix, and investigate its implications. The procedure is the following: 

• First, we perform the computation of all needed quantities, in terms of the trial functions 



B and C. These are the expectation values (3.20) that build up the first cumulant Ai, 
and the normalisation factor mo. 



• Second, we impose the variational equations (3.16) to the trial functions. 

• Third, we read out from those the scattering phases and the T-matrix. 

• Fourth, we obtain the first correction, which is the second cumulant A2, 

• And fifth, we have a look at the high-energy behavior of our results. 

First Cumulant. The first cumulant Ai is the first central moment {S — St). It leads to the 
following approximation 



m J 



(4.1) 



with 

"^^(b, K) = F[St] - 1 = mo exp z (5 - St) - 1 

= moexp z(xK)exp(-2((B,v) + (C,w))) - 1. (4.2) 

The quantity mo is the normalisation factor 

j DvPw e*^' = j PvPw exp [i^ [(v, v) - (w, w)] + i (B, v) + i (C, w; 

= exp('-^[(B,B)-(C,C)]] (4.3) 



Let us consider first the term in (4.2) that involves the functions B and C. This expectation 
value can be written as 



((B, v) + (C, w)) = — / PvPw [(B, v) + (C, w)] exp [iSt{^, w)] 
j PvDw exp y y (v' - w^) + la [(B, v) + (C, w)]^ 



—i d 
rriQ da 



This last expression is simply the normalisation constant mg, with the trial functions multiplied 
by the additional factor a. One obtains 



((B,v) + (C,w)) = -—^ 

niQ da 



^exp(^-— a^[(B,B)-(C,C)]J (4.4) 

-i[(B,B)-(C,C)]. (4.5) 

m 
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Therefore, formula (4.2) for the first term in the cumulant expansion becomes 



Tt'\h,K)=exp [t (xK)]exp 
The first cumulant is thus given by 



2m 



[(B,B)-(C,C)] -1 



Ai = (xk) + ^[(B,B)-(C,C)] 



(4.6) 



(4.7) 



The variational equations are now easily read out. Since we want our expression (4.6) to be 
stationary with respect to the trial functions, it must hold 



-m 



respectively 



C(t) = m 



H xk) (B,C) 
6B{t) 

6C{t) 



(4.9) 



Computation of (xk)- This is done most easily by Fourier expanding the potential function, 
(Xk) = -^J ^v^w / (^^(P) / e^PKi^+'^^W-^-W]. (4.10) 

Since and are linear functions of v, w respectively, the functional integrations can indeed 
be exactly done. Let us first have a look at the v integration : 

J V\ exp ^^"^y (v, v) + i (B, v) + i^p J ds v(s)sgn(t - s] 



lim exp 

T-s>oo 



I 

2m 



+T 



ds 



-T 



1 2^ 



B(s) + -psgn(t - s) 



(4.11) 



The w integration is similar and results in 



lim exp 



T-s>oo 



2m 



ds 



C(s) 



1 



-psgn(-sj 



(4.12) 



From the terms in the squared brackets, the B^ and terms cancel against the normalisation 
constant. Both terms that are seemingly infinite cancel each other, and we are left with 



(Xk) 



dtV {h+—t — 

m 2m 



ds B(s)sgn(t — s) + C(s)sgn(— s) 



From now on we will write x(t) for the trajectory 

J B(s)sgn(t — s) + C(s)sgn(— s) 



, K 1 

h + —t 

m 2m 



x(t). 



Our result is thus 



(Xk) = - j dtVi^it)). 



(4.13) 
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4.1. Variational Equations and the Variational Trajectory 

We see now clearly that the role of the trial functions B and C is to determine a better 
reference path that the particle travels along. As expected, setting B and C to zero reduces to 
the straight line trajectory, and subsequently to the eikonal approximation, here in the variant 
of Abarbanel and Itzykson (AI) |T2], where K appears everywhere in place of k. We will see 



that in our case too, where the choice of trajectory follows from the variational equations (4.8) 



and (4.9), the trajectory can be readily interpreted. The variational equations translate to 



Bit) = -^jdsVVi^is))sgn{s-t), 



C{t) 



which, with the help of the identity 



ds Vl^(x(s))sgn(-t), 



implies in turn 



ds [sgn(s - t)sgn(s - t') - 1] = 2\t - t'\ 



x(t) = h+—t - — I ds VV(x(s))\t - s\ 
m 2m 



(4.14) 
(4.15) 

(4.16) 
(4.17) 



This identity is none other than a integral representation of a solution x(t) to the Newtonian 
equation of motion for a point particle. Indeed, the first two terms on the right hand side are 
the most general free solution, and the third term is the convolution of the gradient of the 
potential function with the Green function of that equation of motion. This is formally seen 
by applying twice a time derivative on both sides of the equation. Since 



It 



\t — s\ = sgn(t — s) and ^sgn(t — s) = 25{t — s). 



(4.18) 



the computation is straightforward : 



mi(i) = -V\/(x(t)). 



(4.19) 



Interestingly enough, in this approximation of quantum mechanical scattering theory, an im- 
portant feature of which is the scattering integral equation, its classical counterpart appears. 



Asymptotic Behavior. In order to get a feeling for the role played by the two fixed pa- 
rameters b and K in this trajectory, let us consider more closely its asymptotic behavior. We 
will always assume in the subsequent considerations that the potential falls off rapidly enough 
so that the integrals involved converge. For very large positive times t, the term t — s m the 
integral can be considered always positive. One can then rewrite the trajectory as 



x(t) *^ b/ + tv/. 



(4.20) 
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where the starting point bj and asymptotic final velocity vj are related to b and K by 

hf := h + ^ [ dsVV{x{s))s :=h + lh (4.21) 
2m J 2 

Vf := -- — f dsVV(^(s)) :=- + -^. (4.22) 

^ m 2m J ^ ^ m 2m ^ ^ 

In the very same vein, one obtains for very large negative times 

x(t) b, + tvi, (4.23) 



with 



b - ^b (4.24) 

^ _ _9_ 

m 2m 



(4.25) 



These relations implies immediately 



b = ^(b,, + b^) (4.26) 
^ J(v, + v^). (4.27) 



m 

One deduces from these relations that the integral equation for x(t) is (at least formally) 
equivalent to the boundary value problem given by the Newtonian differential equation of 
motion, with the boundary conditions determined through 

b = - lim {x(t) + x(-t) - t [i(t) - i(-t)] } (4.28) 



Equation (4.27) says that K keeps very nicely its meaning of mean momentum, something 
that had a priori no reason to be expected from the variational principle. The quantity Q, 
on the other hand, is then the momentum transfer due to the defiection in this trajectory 
only, but cannot be compared with the momentum transfer q of the scattering process we are 
investigating. Since the energy is conserved along the trajectory, Q is a purely transversal 
quantity. The interpretation of b for a general potential is on the contrary more difficult. Sure 
is that if b is large enough, the trajectory is nothing more than the straight eikonal line. This 
follows from the fact that if the potential at b and beyond is negligible, then this straight 
trajectory is a solution to the integral equation. Nevertheless, since one integrates over b to 
obtain the scattering amplitude, we see that every classical trajectory with mean momentum 
K contributes in our approximation to the scattering process. It is very similar in spirit with 
the eikonal expansion, where all straight line trajectories with velocity k/m are taken into 
account. 
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Constants of Motion. The asymptotic conditions we derived above permit us to derive 
simple expressions for the energy and, in the case of a potential with rotational symmetry, for 



the angular momentum of this classical curve. According to (4.23) and (4.20), the angular 
momentum mx A x at infinity is given by 



L = mhf A V/ = mhi A Vj. 
It is possible to combine these two to obtain 



m 

~2 
m 



{hi A Vi + b/ A V/) 
K 



, , b - -b ) A , 

2 VV 2 y \m 2m 

b A K + -b A Q 

4 



Q 



Ih 1 A 



The same procedure applied to the energy gives as a result 

E 



K2 Q2 

h — . 

2m 8m ' 



K Q 

m 2m 



(4.30) 

(4.31) 
(4.32) 
(4.33) 

(4.34) 



These two relations prove to be useful for instance if one is interested in a high K expansion. 



Symmetries. Under some assumptions on the potential, the trajectory exhibit some further 
symmetries. Let us suppose that the potential function is such that it depends only on the 
modulus of the components that are perpendicular and parallel to the mean momentum. 



y(x) 



F(xJ,xl 



(4.35) 



To this class of functions belong of course the potentials that possess rotational symmetry, 
but not only those, since each component could be treated differently, like in a " cylindrically" 
symmetric potential. Although this last category may look a bit artificial, it has some im- 
portance, since we will encounter such a function in the ray representation of the T-matrix, 
where the perpendicular and longitudinal parts are handled in different ways. Anyway, for 
such potentials it clearly holds that 



VxV(x) = /(xJ,xi)xx and V,, V(x) = (^(xj, xi)x 



(4.36) 



for some scalar valued functions / and g that may or may not differ. We claim that this 
implies in turn the symmetry properties 



x^(— t) = x_L(t) and X||(— t) = — X||(t). 



(4.37) 



These properties, together with the conditions ( 4.36[ ) we stated above, are clearly seen to be 
compatible with the defining equation (4.17). Although physically intuitive, to affirm that 



x(t) must possess these, one can argue as follows: the standard way to approach the solution 
to such an integral equation is to proceed through iteration. Starting with 



xo(t) 



b + — t, 
m 



(4.38) 
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one defines recursively 



x„ t := b + -t - — 
m 2K 



j dsVV{^n^i{s))\t 



(4.39) 



with the velocity v = K/m. In a high energy hmit, i.e. fixed velocity but high K, the so- 
defined x„(t) are the successive approximations to x(t) up to order K~'^. Since ^oit) possesses 
the symmetry properties stated above, by induction over n all terms satisfy them. Since such 
symmetry properties cannot depend on the very magnitude of K, x(t) must in general obey 
these, too. 

This has as a consequence, for example, as can be easily checked, that x(0) is purely perpen- 
dicular, and x(0) purely longitudinal, implying, if we denote |x(t)| with r(t). 



d 



r(t) 



1 



t=o 



r(0) 



x(0) ■i(O) = 0. 



(4.40) 



In other words, x(0) is the closest approach to the scattering center. 



4.2. The Scattering Phases 

We have computed in the last section the expression for the T-matrix in first order in the 



cumulant expansion, in terms of (xk)? in equation (4.6). When B and C do satisfy the 

1 (4.41) 



variational equations (|4.14|) and (|4.15|), it becomes 

exp {iXo + iXi 



rf~^)(b,K) 



where we introduced a bit of notation : 



Xo(b,K) = - j dtV{^{t)) = {xK) 



(4.42) 



and 



Xi(b,K) = -^ J dt J c/s VV(x(t)) ■ VV(x(s))|t-s|. (4.43) 



For potentials with rotational symmetry, since b is perpendicular to K one has 

Tf"^)(b,K) =Tf~^)(6,ir). 



(4.44) 



A first, very rough estimate of Xq for small values of b can be done as follows: we denote the 
scale over which the potential takes on sensible values, with average value Vq, by R; since the 
particle following this trajectory would cross this region in a time ~ R/v, where v is the mean 
velocity K/m, one has roughly 



R 



-e KR, 



Vo/Kv. 



(4.45) 
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In the last equality we have introduced the dimensionless quantity e, which describes the 
relative importance of the potential energy with respect to the kinetic energy. A similar 
argument in the case of Xi shows 

2 



a /AV 



K 



(4.46) 



where a stands for the scale over which the potential varies appreciably, by an amount AV . 
In the case where the quantities a and R are essentially the same, one has 



-e' KR. 



(4.47) 



Of course, a and R can sometimes differ considerably, as in the case of a Woods-Saxon potential, 
that looks like a well with rounded edges. In that particular case, 



-e^ Ka. 



(4.48) 



The larger the energy of the incoming particle, the smaller epsilon, the more dominant is Xq 
with respect to Xi. 

These two phases are very similar in shape to those of the eikonal expansion, up to first order 
in inverse momentum, except apparently for the minus sign in front of Xi. The first two phases 
of the eikonal expansion of Wallace read for a general potential 



dtV \ h + K—t 

m 



+ 



1 

Am 



dt I dsVV [h + K—t ] - VV [h + K—s ] \t - s\. (4.49) 



m 



m 



We will nevertheless show later that they coincide in the high energy limit. 

However, the two scattering phases Xq and Xi satisfy an amusing relation that has no eikonal 

counterpart. It is indeed possible to show that 



Xo - Xi = K ■ (b 



/ 



Kb, 



(4.50) 



where we have used the notation introduced for the asymptotic properties of the trajectory in 
(4.20) and (4.23). The proof reads as follows: Considering the defining equation (4.17) and 
its classical property (4.19), one can first rewrite Xi as 

dtVV{^{t)) ■ 



Xi 



x(t) - b 1 

m 



m 
' 2" 



dt x(t) 



X t) - b 1 

m 



(4.51) 
(4.52) 



Using finite bounds for the integration interval, one can integrate by parts and consider each 
term on its own. 



Xi|t = -^ x(t) 



x(t) - b 1 

m 





t=+T 








m f 




+ 






t=-T 


^ J- 



+^ K 

rftx^(t) x(t) 

2^ m 



t=+T 
t=-T 



The second term on the right is simply 



m 
~2 



[ dt±\t)= [ dt [E -V{x{t))] = 2TE + Xo\j,. 

J-T J-T 



(4.53) 



(4.54) 
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By inserting the asymptotic trajectory into the two other terms, one sees after a small calcu- 



lation using (4.34), that all dependence on T disappears, and what is left is the formula (4.50) 
claimed above. 



4.3. First Correction to the Variational Approximation 

In this section, we want to go a step further and compute the second order term A2 in the 
cumulant expansion of (e*^"^). However, to keep things simple we will stick to the solutions 
of the variational principle we derived earlier. Especially, we expect the appearance of a real 



term. The main results of this rather long and technical section are equation (4.84) combined 



with (4.78), and equation (4.88). 



This second order term is the second central moment, or variance (see (3.26)). It is related to 
the first two moments by 

X2 = (AS^) - {ASf. 
This means that our approximation of the integrand in the T-matrix in is going to be 



(4.55) 



p(3-3) 



.K 



m 



J c/26e-'^-'iTf-^^(b,K), 



with 



.(3-3) 



(b,K) 



moexp ( iXi 



A2 - 1. 



As usual, tuq stands for the normalization factor 



mo 



exp 



--|(B,B)-(C.C)| 



(4.56) 
(4.57) 

(4.58) 
(4.59) 



To this aim, the six terms entering 

{AS') = ((xK-(B,v)-(C,w))^> 

= (x^> + ((B,v)^> + ((C,w)^> 

- 2 (xk (B, v)) - 2 ixK (C, w)) + 2 ((B, v) (C, w)) 

need to be calculated. The last term causes no worries : since v and w are uncorrelated, the 
expectation value factorize out, and the term will cancel its corresponding part in (AS)'^ in 
A2. 

The four terms involving the trial functions B and C can be computed using a similar procedure 
we did when we calculated the first cumulant: we replace B, or C respectively, with the help 
of a fictitious variable a in the action St by aB, or aC, and derivation of mo with respect 
to a brings down each time a factor i(B,v). More concretely, we have for instance, for each 
positive integer n, 
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((B,vr) = — 

mo 

1 d"" 



tuq \i J da"- 
1 rf" 



VvVw (B,v)" e^?^'^*(^'-^')+'^(^'^)+*(*^'^) 



mo V* 



a=l 



a=l 



exp<j-— [a2(B,B)-(C,C 



(4.60) 
(4.61) 
(4.62) 



An even more concise form can be written after recognizing that this very last expression 
represents a rescaled Hermite polynomial Indeed, we obtain 



((B,v)") = z" 



fB,B 



m 



n/2 



Her 



m 



where Hen{x) is the Hermite polynomial of degree n, defined [TT] by 

HeJx) = (-ire^"/'^e-^"/^ 



dx^- 



In a completely analogous way, one gets 



((C,w)'^) = z" (C,C 

m 



n/2 



Her 



m 



■(C,C) 



Since He2{x) = — 1, it follows immediately that 



;(B,v)^> 



1 



m 
i 

-( 
m 



(B,B) + — (B,B) 



m^ 



B,B) + ((B,v))' 



respectively 



((C,w)^> 



■-(C,C) + — (C,C)2 
m m^ 

' (c,c) + ((c,w))^ 



m 



In a similar vein. 



(Xk(B,v)) 



— / V^rVw XKiB,^r)e'^' 
mo J 



1 d 



mo (id 
, 1 d 
rriQ da 



a=l 



a=l 



(xk) ((B,v))-^ 



d 
da 



(4.63) 



(4.64) 



(4.65) 



(4.66) 
(4.67) 

(4.68) 
(4.69) 

(4.70) 
(4.71) 
(4.72) 
(4.73) 



a=l 
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In the last three hnes, the upper script " denotes the usual quantity with B replaced by aB. 
We know already what (xk) is from our study of the first cumulant: 



(Xk)'' = - / dtV[h + -t-^ 

m 2m 



ds [aB{s) sgn(t — s) + C(s) sgn(— s)] 



It follows then 

iXK (B,v)) = (xk) ((B,v))- 
and in an analogous way, 

(Xk(C,w)) = (xk) ((C,w)) 



2m 



J dt j ds VV(x(t)) ■ B(s)sgn(t - s), (4.74) 



2m 



dt / dsW{x{t)) ■ C(s)sgn(-s). 



We now gather equations (|4.67|), (|4.69|), (|4.74|), (|4.75|) to obtain 

A2 = 



(x^>-(Xk)V-(B,B) 



feci 



m 



m 



m 



+ — dt rfs V\^(x(t)) ■ [C(s)sgn(-s) + B(s)sgn(t- s)] 



(4.75) 



(4.76) 



(4.77) 



If B and C are taken to satisfy the variational equations we derived earlier, we obtain the 
representation 



A, 



[x'k) - {XKf + VV{^{t))VV{^{s)) \t - 



(4.78) 



As expected, divergent phases again cancel with the help of the antivelocity. Therefore, it 
holds 



Tf-^)(b,K) 



mo exp ( i\i - -As 



(4.79) 



= exp(^zXo + 2zXi-^[(x^>-X2]j -1 
We are thus left with the task to compute (xk)- 

Computation of (xk)- Again, this can be can be done by Fourier expanding the potential, 
and solving the v,w functional integrals : 



1 

mo 
1 



dti j dt2 

2 



- / F ( b + — t + x^(t) - x^(0) 
m 



(4.80) 



mo 

■ exp <{ z ^ p 



(ipi [ dp 



i=l 



(27r)3 J (27r)3 
K 



PvPw e^^* 1^(pi)l^(p2 



b -|- tj -|- x^(tj) 'x.^iti 
m 



(4.81) 
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Let us consider first the v integration. 



/ 



"Dvexp 



exp 



2m 



-(v • v) + i (B, v) + - Pi / ds v{s)sgn{ti - s) 

i=i 

J ds (b{s) + ^ 5^ Pi sgn(ti - s) 



(4.82) 



Similarly, the w integration results in 



exp 



2m 



J ds ^C(s) + ^^p, sgn(-s)j 



(4.83) 



We see once again that the antivelocity cancels the divergent phases. From the squared terms, 
and are killed by the normalisation constant, and the others combine to give the result 



(Xk> = dt, dt 



dpi f dp 



^(Pl)^(P2) 



(27r)3y (27r)3 



• exp <! i [pi • x(ii) + p2 • x(i2)] exp ( — pi • pa \ti - ^2! 



(4.84) 



It is also possible to integrate over the momenta to obtain a real space representation. To this 
aim, we first rewrite it as 



<X^> = JdtiJ dt, J ^ V{p)V (^x(tO + ±^p\t, - t2|) e 



ip-x(t2) 



(4.85) 



and expand the potential function in a power series. 



It suffices then to note that p"y(p) = (-i)l"lL>"y(p) to conclude 



\a\>0 



\a\ 



(^K> = E ^ j / / L'"y(x(iO)^"^(x(i2)) 1^1 - t2|l"l. 



(4.87) 



Obviously, the first term in this infinite sum is (xk)^, and those with \a\ — 1 build twice 
the phase Xi. Therefore, the second cumulant can also be written as 

A2 = E ^ (^) [dti [ dt2 L'"V^(x(ii))L>"\/(x(i2)) 1^1 - t2|l"l, (4.88) 



|a|>2 



where the sum starts at |q;| = 2 only. 

While this form is useful for a high energy expansion, it is, on the other hand, rather unpractical 
for numerical computations. 
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4.4. High Energy Expansion 

(3—3) 

It is interesting to investigate the large energy behavior of T-;_^j^ , so that it can be then im- 
mediately compared to the systematic expansion done by Wallace in |1] . In such an approach, 
all quantities are expanded in inverse powers of k, while the velocity term k/m is treated as 
fixed, as appropriate for a high energy situation. Starting from the phases Xq and Xi and 
the second cumulant we just derived, we will investigate the first and second order of the high 
energy expansion, and obtain the corresponding T-Matrices, defined through 

Tt'\h) = exp (z [xo(b) + xi(b)]) - 1, (4.89) 

and 

Tl'j-'\h) = exp (z [xo(b) + xi(b) + X2(b)] - uM) - 1. (4.90) 

The subscripts under the different terms indicates the power of 1/A; attached to it in the 1/k 
expansion. In doing so, one has to be careful with the corrections arising from the fact that it 
is not the initial momentum k, but the mean momentum K, that appears in our path-integral 
formulation of the T-matrix. However, since 



these corrections will not appear at first order. This is why notably no imaginary phase uJi{h) 
can appear at all at first order . We will see later how to handle these corrections. 

First Order. The leading term of the trajectory x(t) in a high K expansion is given by 

Xo(t) := b + —t. (4.92) 

m 

The first order contribution is obtained in the following way: 

K If 

x(t) = h+—t dsW(:>c(s))\t - s\ 

m 2m J 

= xo(t) j dsVV (xo(s)) \t-s\ + --- , (4.93) 

where we wrote v for K/m. This implies in turn that the leading phase Xq can be written as 
Xo = - j dt Vi^it)) = - j dtV (xo(t)) + 

+ ^ j ds j dtVV (xo(t)) ■ VV (xo(s)) \t-s\ + --- . (4.94) 



It is seen that the phase Xi to first order is simply given by 
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ds / dt VV (xo(t)) ■ VV (xo(s)) \t - s\ + 



(4.95) 



and that the second cumulant contains only terms of order and higher. As we have just 

argued, the difference between K and k plays here no role. Therefore, the leading phase Xo(b) 
and the first order phase Xi(b) in the high-/c expansion are given by 



Xo(b) 



dtV [h 



and 



K 



m 



(4.96) 



XM = 1^ jdt I dsVVih+^t).Wv(h+^s]\t-s\. 



For spherically symmetric potentials, Xo(b) can be immediately rewritten as 

1 



Xoib) 



dZ V (r) , 



while we show in Appendix B.l that 



(4.97) 



(4.98) 



(4.99) 



In both of these expressions, the argument r of the potential is meant to be \/W^r~Z^. These 
two phases coincide exactly with the first two terms of the eikonal expansion from Wallace. 
Thus the first cumulant of the variational procedure contains both the leading order as well 
as the first order of the high energy expansion. As we will see, it contains even the imaginary 
part of the second order as well! 



Second Order. The second order term is a bit more subtle. There are three sorts of con- 
tributions to it. 

First, the two cumulants Ai and A2 contain terms of that order. While Ai contains only a 
complex part, that we will denote X2^~^''(b), the second cumulant contains to that order only 



a real part, which we call uj^ ^^b)- Indeed, the expression (4.88) we derived for the second 
cumulant can be rewritten as 



E ^ (^)' '/ / dt,D-V{^{t,))D-V{^{t,)) |ti-t2|l"l, (4.100) 



« ^ I / 1,1 1 \ 

\2 = 

\a\>2 



with V = K/m. The first term of this sum, with the straight line trajectory Xo(t) as argument 
in the potential functions, is already second order in and thus also in l/fc, and is purely 
real. All other terms are of higher order. 

Second, the fact that the factor iK/m stands in front of the T-matrix, and not ik/m, gives 
rise to a second order contribution. We will denote the imaginary and real contributions with 
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(b) and (b) respectively. 
Third, the term of order zero in our eikonal representation was found to be (see equation 



(4.94)), after the substitution Z = tK/m, 

-'I^JdZV{h + ZK) = -j(^l-^ + --^ ' jdZV{h + ZK). (4.101) 

Although we could ignore the difference in first order, it must be taken into account in second 
order. We will call the imaginary contribution u}2^{h) and the real one uj2^{h). 

Computation of x^2~^'^ • definition this is done by expanding the complex quantities 
Xo(b) and -^i(b) up to second order. Let us first consider Xi. Since Xiih) is already first 
order, we need to expand once around the leading trajectory xo(t). 

Xi = dsW{x{t)) ■W{x{s))\t~ s\ 

1 
4m 



^jdtjds Vy(xo(t)) ■ Vy(xo(s))|t-s| 
-2^ / duViVjV{xo{t))VjV{xo{s))ViV{xo{u))\t - s\\t - u\. 
Similarly, one has to expand Xq to get the second order contribution. Starting from 

Xo = - [ dtV{x{t)), (4.102) 



one gets, after expanding around xo(t) to second order, 

Xo^- j dt F(xo(t)) + ^jdt j ds VV{^o{t)) ■ VV{^{s))\t - s\ 

~\i^2rri) j ^^^^^"^»'^i^(^o(i))Vil^(x(s))VjK(x(u))|t-s||t-M|. 

The second order term is then given by the third term in the last equality with xo(s) and xo(m) 
instead of x(s) and x(m), while the second term must be again expanded around xo(s) and 
gives rise to another contribution of the very same form. After getting the numerical factors 
right, the total contribution from Xq and Xi is 



^(3-3) _ 

X2 g^2 



j dtdsduViVjV{^Q{t))VjV{i^Q{s))ViV{^Q{u))\t-s\\t-ul (4.103) 



with V = K/m. For spherically symmetric potentials, it is shown in the appendix B.3 that it 
reduces to the form 

(4.104) 



1 

8^2 
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Computation of ^nd Wg^. These two terms are the complex and real contributions 
coming from the K factor in front of the T-matrix: 



T, 



(3-3) 



.K 
m 
. k 



A, 



d% l- 



e "^"^^ I exp 



- 1 



m J \ Sk"^ 

To absorb this term into a correction phase depending on b only, we observe that 



q^ 



1 - ^ 1 exp (-iq . b) = + ^ J exp (-zq ■ b) , (4.105) 
with Ah being the 2-dimensional Laplacian operator with respect to b. It is now possible to 



integrate by parts and to incorporate the correction into the phase. In the appendix D. 
show that the final results are 



we 



and 



^2Hb) = — (V.xo)^. 



(4.106) 
(4.107) 



Computation of and uj2^ . As explained above, these two terms appears because of the 
factor of K present in the leading order phase in the eikonal representation. As for and 
uj2^ , they come with a factor that can be transformed away by a Laplacian operator acting 
on the scattering phases after an integration by parts. Detailed calculations are provided in 
appendix inn The results are 



X^^(b) 



8A;2 



((V6Xo)\o - AbXo) , 



and 



c^2^(b) 



8P 



(2 {VbXof + XoAfeXo) 



(4.108) 



(4.109) 



Computation of ui2^ We already mentioned that 1^2^ ^' is the first term in the infinite 



(3-3) 



sum (4.88) that describes the second cumulant. Getting all the numerical factors right, it thus 
reads 



(3-3) 



(b) 



16i^2 



J dt J ds\t- s|2 ViVjV (xo(t)) ViVjV (xo(s)) (4.110) 



Appendix B2 provides the derivation of its reduced form for a spherically symmetric potential. 
There it is shown that 



(3-3) 



1 



8K^ 



2Xo^ + Xo [Xo + bXo] + -^XoXo 



1 



(4.111) 
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It suffices now to sum over all the contributions to obtain our second order scattering phases. 
We have, written here explicitly only for potentials with spherical symmetry, 



b 



with AbXo = x'o + x'o/b. For the imaginary term, 

X2ib) = xt'\b) + x^'ib) + x'2'ib) 



5 d 
3 db 



1 

-b^ — 
3 db^ 



'3 



dZ V^(r) - 



bX6 
24A;2' 



(4.112) 



(4.113) 



These expressions turn out to be exactly the same as those in the expansion from Wallace (see 



expressions (2.45) to (2.48)). It is remarkable that the variational procedure, in the very ffist 
cumulant, contains already all imaginary terms up to the second order, plus a part from the 
real term uj2- The ffist correction to the variational procedure, i.e. the second cumulant, then 
completes it. 



4.5. Unitarity of the Approximation 

The investigation of the unitarity of the scattering matrix, or the validity of the optical theo- 
rem, in our variational approximation is not an easy task. Even at the very beginning, although 



(2.32) and (2.41) are exact representations, its validity is but self-evident. We have seen that 
the variational approximation is close to an impact parameter representation of the T-matrix. 
Especially, we just showed that it reduces to the eikonal expansion for high energies. We will 
try to gain some insights in this question within a very brief look at the unitary of impact 
parameters representations. We will consider in the following central potentials only. 
We denote by x the imaginary and by u the real part of the scattering phases, and by S(b) 
the quantity 

S{b) := e^'^^'')-'^^''), (4.114) 
so that the T-matrix elements in the impact parameter (eikonal) formalism are given by (see 



equations (4.89) and (4.90)) 



T,^/ = i- [ d% e-'""-^ T(6), T{b) = S{b) - 1. (4.115) 

171 J 

The optical theorem reads 

^tot = y»/(0)- (4.116) 

It relates the total cross section to the imaginary part of the scattering amplitude /, and 
follows directly from the condition of unitarity of the scattering matrix (a derivation can be 
found, for example, in pQ). The scattering amplitude is in turn related to the T-matrix through 

rn 

m) = --T.^f. {h=l) (4.117) 
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When T is written in the form (4.115), the following relation must therefore hold: 

-2 / (fb [^(S) - 1] . 



<ytot 



The total cross section is, on the other hand, by definition 



Ctot 



(27r)2 

Since q = 2ksm{9/2), it holds that 



d% I d%' e-''^^e''^^'T{b)T*{b' 



qdq = 2k^ sin - cos - = sm{e)d9 



with implies in turn 



dn{e,, 



d^q. 

q<2k 



(4.118) 

(4.119) 
(4.120) 

(4.121) 
(4.122) 



The range of integration is restricted to the physically possible momentum transfers whose 
value does not exceed 2k. However, one can also integrate over all possible values of q, and 
then subtract the non-physical scattering cross section, for which q > 2k: 



C^tot 



(27r 



^ J d\J d% J d%' e-''i-'^e*'i-^'T(6)T*(6') " ^np = ^ d% \T{b)f - a^p, (4.123) 



where cr„p is 



cr 



np •- 



(2vr) 



d\ 



q>2k 



d% e-'^'-^Tib) 



> 



If the optical theorem is to hold, the following relation must thus be fulfilled: 

d% {\T{b)f + 2mS{b))-l]} 

'd% [\S{b)\'-l] 



np 



d% [t 



-2uj{b) 



(4.124) 

(4.125) 
(4.126) 
(4.127) 



Since Cnp is positive definite, the real part uj of the scattering phases is thus seen to play an 
essential role in the theory. At high energy, we have seen that the real term uj predicted by 
the variational approximation is the term(2.48), 

6^2 (ft) 



iPXo) ■ 



16A;2 b db 

The non-physical cross section can then be calculated as 



np 



-> -2 



d%UJ2{b) 



An 



16k' 



Jo '^^db^^'^' 



(4.128) 

(4.129) 
(4.130) 
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This expression vanishes for potentials that are regular at the origin, but this contribution is 
finite, for instance, for a Yukawa potential. The question how well unitarity is fulfilled at high 

energy thus depends strongly on the specific potential. 

One of the main differences between the variational approximation in the eikonal representation 
and in the ray representation is, as we will see, that in the latter, such a real part is contained 
in the approximation, while we have seen that in the former it appears in the second cumulant 
only. It should therefore not be considered as a surprise that the variational approximation in 
the ray representation will turn out to be superior than its eikonal counterpart. 
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5. Variational Approximation in the Ray Representa- 
tion 

Besides the facts that the ray representation does not treat the parallel and perpendicular 
components on an equal footing, and that the trial function C{t) is yet only a scalar, all 
calculations in the ray representation resemble very much those we have performed in the 
eikonal representation. However, we will encounter here new interesting features, that render 
this approach very different from the eikonal one of the last section. We will proceed in the 
very same order, to obtain first, the approximation 

T^^-^^ =1-1 d^h e-^^"i "')(b, q, K), (5.1) 
m J 

that contains only the first cumulant, and, second, 

T/i/ =i^j ^i26e-^-''Tf-^)(b,q,K), (5.2) 
with the first correction (second cumulant) included. 



First Cumulant. Obviously, the first cumulant is given by the suitably adapted formula 

(5.3) 



(ra, i.e. 



Af-'^ = (W + ^[(B,B)-(C,C)], 



and the variational equations for B and C are the equivalent of (4.8) and (4.9): 

5{Xray) (B, C) 



B(t) = -m 



respectively 



C(t) = m 



6{Xray) (B, C) 

6C{t) 



(5.4) 



(5.5) 



Computation of (Xray)- Since the computations are very similar to those we already have 



done, they are only reported in appendix A. There it is shown that the expectation value 
(Xray) has the following form: 



iXraylay = - I dt I J^,V (p) CXp ( ^p (b + U exp ( - — |t|pi 



'^^^ (~2m^ {/ ^^^^ ^^^^^^ - ^) - sgiul-'S)] + C(s)sgn||(-s)|^ 



(5.6) 



The symbols _L (or ||), such as in sgn_|_(s), have the meaning that it multiplies only the 
components of the 3-dimensional vector that is involved in the product that are perpendicular 
(or respectively parallel) to K (for example psgn_|_(s) = p_LSgn(s)). We will use this convention 
in the following, as it permits to render most expressions more concisely and elegantly. 
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This expression for (Xray) contains something new with respect to the (3-3) version. The 
exponential factor that contains a quadratic term in does not permit to return immediately 
to real space, as we could do in the eikonal representation. Nevertheless, we interpret the 
following quantity K-(j) as being a new, effective potential, defined in Fourier space through 
the (Fresnel) Gaussian transformation 



V.it){p) ■■=V{p) exp(-^a(t)pl 



with a(t) = — \t\. 

m 



(5.7) 



This new potential, in real space, is the convolution of the potential V with a Gaussian filtering 
function with typical filter size ^/o^. However, the filter being complex, this effective potential 
is a priori also a complex quantity. This is an essential feature of the ray representation. Note 
that at time t = 0, the effective potential is the same as the potential V, and that is also 
vanishes in the limit of asymptotic times. 

The arguments that multiply the linear term in p can be then defined as the trajectory 



= b + [ds B{s) [sgn(t -s) + sgn^(s)] - KC(s)sgn(s) 

m 2m [_ J 

and the final result for (Xray) has a similar form as in the eikonal representation 

iXray) = - J dt V^(t){z{t)) . 

5.1. Variational Equations and the Variational Trajectory 



(5.^ 



(5.9) 



We can now proceed to solve the variational equations in the ray representation. Just as in the 
(3-3) case, setting B and C to zero reduces the trajectory to the straight lines h + tprayit)/m. 



However, if these are taken to satisfy the variational equations (5.4) and (5.5), it must first 
hold that 



Bit) = -ljds VK(s)(z(s)) [sgn(s - t) - sgn||(-t)] 



and 



C(t) 



^Jds VK(,)(z(s))sgn||(-t). 



and, second, after inserting these expressions into the trajectory z, we get 



z{t) =h + -t+ ^\t\ ds VK(s)(z(s)) \\t - s| - |t| 

m 2m 2m J 



(5.10) 
(5.11) 

(5.12) 



This trajectory is a complex trajectory, and possesses some properties analogous to those in 
the eikonal representation. The component of the trajectory that is parallel to K also behaves 
classically. We have, by differentiating twice. 



z||(t) 



1 



m 



V||Kw(z(t)). 



(5.13) 
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However, the perpendicular part gets a kick at t = 0. Indeed, the two terms that contain \t\ 
give rise to 26(t), and we get 



i±{t) = --VxK{t)(z(t)) + (- + - I rfs VxKw(z(s)) 
m \mm' 



(5.14) 



Asymptotic Behavior, Constants of Motion: Symmetric potentials. The effective 
potential \4-(t)(z) has a special dependence on the perpendicular components of z. However, 
it still depends only on the modulus of its parallel and perpendicular components, if V does 
so, although in different ways. In that case, by the very same argument as in the eikonal 
representation, the symmetry properties 



Z±(-t) = Zx(t), Z||(-t) = -Z,|(t) 



(5.15) 



still hold, although Va(t) is generally not spherically symmetric. The defining equation (5.12) 
can be then rewritten such that the unnatural subscripts are eliminated. 



z{t) = h+-t + ^\t\-— ds VK(.)(z(s)) [\t -s\- \t\ 
m 2m 2m J 



(5.16) 



The trajectory (5.16) exhibits interesting properties. First, for very large times, it holds that 

I f— >itoo 



t-s 



t 



(5.17) 



so that the integral part of the trajectory goes to a constant. We name its parallel component 
Bray/2, and its perpendicular Cray, 



B 



ray 



''ray 



SO that 



Z(t) *^ b + Cray ± ^B, 

This implies immediately 



m 
1 
2m 

1. 



ds VK{s)(z(s))s 
ds VK(s)(z(s))|s| 



_q^, ray 

™^^2"™^ m 2m ^'^ +^^/.*- 



. .... i->±oo K q 

zvtj — > — ± 

m 2m 



kf/m t — )• +00 
kj/m if: — )■ — oo 



(5.18) 
(5.19) 

(5.20) 
(5.21) 



The asymptotic momenta of the variational trajectory (5.16) are thus precisely those of the 
scattering process, kj and kj. Both the mean momentum K and the momentum transfer q of 
the scattering process have now the same meaning in the variational trajectory. 
Furthermore, the "energy", which is now a complex number. 



E{t) = ]^mi\t) + V^(tMt)) 



(5.22) 
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is at infinity 

E{±oo) = ^ = ^, (5.23) 
2m 2m 

since the potential vanishes. This is precisely the scattering energy of the process we are 
investigating. However, there is no conservation of energy for finite times, since the effective 
potential has an explicit time dependence. Nevertheless, the kick at time t = is elastic. This 
can be seen from the fact at t = the perpendicular components of z{t) get reversed, while 



the longitudinal part stays untouched: from (5.16) and (5.15) it follows 



= ^ ± ^ - i/^^ V,y„„=g„(-,) r^^jds V,14,.„ (5.24) 

which is to say 



Also, it holds 



Z||(0+) = Z||(0_), and i±{0+) = -i±{0^). (5.25) 



z(0) = b. (5.26) 



A particle following this trajectory thus starts with the right energy of the scattering process, 
enters the interaction region where its energy changes, and becomes complex (see figure [l] for 
an illustration), at time t = experiences an elastic reflection on the plane perpendicular to 
K, and then gets out with the same energy as when it entered. 




Figure 1: The real (red line) and imaginary part (blue line) of the "energy" of the variational 
trajectory in the ray representation, as a function of time, in units of the scattering energy 
for a particular choice of parameters and a Gaussian potential. With the parameters chosen, 
the target is about 2 fm large, and the incoming velocity about 0.5c. This implies a typical 
interaction time scale, expressed in terms of length, of 4 fm. 
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The analysis of the angular momentum for spherically symmetric potentials also provides 
interesting insights. Since the effective potential itself does not possess rotational symmetry, 
the conservation of the three components of the angular momentum is not achieved. However, 
since V^(t) is still symmetrical in the perpendicular plane, the conservation of its longitudinal 
component is ensured, except at t = 0, where the kick takes place, when it suffers a change. 
This is seen from its definition 

L|| = mzx(t) A zx(t), (5.27) 

(5.28) 



together with the fact (see conditions 4.36) that 

VxKw(z(t)) = /(zl(t),zJ(t))z^(t), 
for a scalar function /, from which follows, for any non-zero time t, 
d 



dt 



L||(t) = mz±{t) A z±{t) = -f z±{t) A z±{t) =0 (t ^ 0). 



(5.29) 



The other components of L would mix in this last equality the functions of / and g of (4.36), 
so that they are in general not conserved. Let us have a look at the change in Ly at time 
t = 0: From (5.25) and (5.26) it follows at once 

Lj|(t = 0±) = mzx(O) A zx(0±) 



±m b A z_l(Oh 



(5.30) 



so that L|| just flips sign at 0. The asymptotic angular momentum can be calculated from 
KM and K27\. It reads 



L||(±oo) = m (b + C 

ray 



q 

2m 



(5.31) 



While the asymptotic energy is the same as in the reference straight line trajectory in the ray 



representation, h + 1 



It^ray (0 



and is thus unaffected by the variational procedure, the angular 



momentum gets a correction in terms of the quantity C 



ray 



5.2. The Scattering Phases 



The scattering phases in the ray representation are obtained from (5.3). We already computed 



the first phase {Xray) with (5.9) as a result. The second phase is found by inserting the 
variational equations (5.10) and (5.11) into the second part of (5.3). It becomes 



Ti^-^\h, q, K) = exp (tX^^''^ + tX^ 



(3-1) 



- 1 



with the definitions: 



Xt'\h,q,K) = - j dtV^^tMt)) = iXray) 



and 



xf~')(b,q,K) = -— / dt I rfs VK(i)(z(t))VK(,)(z(s)) [\t-s\ - |tU - \s\ 



(5.32) 



(5.33) 



(5.34) 
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Again, if the potential V satifies the conditions (4.36), the subscripts ± can be left out. The 
main difference between these two terms and their eikonal counterparts is that they are both 
complex quantities, and contain thus much more information in comparison. In the case of a 
spherically symmetric potential, it holds 



Tt'\h,ci,K) = Tr"{b,q,K,h-q^ 



(3-1), 



(5.35) 



since q and b are perpendicular to K. The relation (4.50) between these two phases has in 
this case the form: 



X, 



(3-1) 



(3-1) 



K ■ 'Bray + Q " Cray + dt 



2m 



k 



Ik, I. 



(5.36) 



The integral on the right is well defined, since, as we have shown above, the integrand takes 
on non-zero values only near the scattering region (see equation (5.23), or figure [T]). The 
derivation is in all aspects similar. After inspection of (5.12), (5.13) and (5.14), the second 
phase can namely written as 



(3-1) 



+ - 



m 
T 
m 



dt'i{t) 

- + - I ds V±K(,) 
m m J 



K 



q 



m 2m 



/ X , K q , 

z(t) - b 1 - —\ 

m 2m 



[z(0) - b] 



(5.37) 



The last term is zero, so that we get 



(3-1) 



m 

T 



dtz{t) 



/ X , K q 

z(t) -b 1 - — 

m 2m 



t\ 



(5.38) 



We use the same trick as we already did to decompose the integral in different terms, 

t=+T 



XI 



(3-1) 



m 



z(t) 



/ X , K q , , 

z(t) - b 1- —\t\ 

m 2m 



m 



K 



+ — / dtz\t) - —z{t) 



t=+T 



t=-T 



t=-T 
+T 



dt z(t)sgn(t). 



(5.39) 
(5.40) 



and recover X, 



(3-1) 



smce 

m 
~2 



+T 



dtz'^{t) 



+T 

■T 
+T 



dt [i?(t)-K(t)(z(t))] 



dt E{t) + X^^~^^ 



(5.41) 



After inserting the asymptotic behaviour (5.20) of the trajectory in the other terms and some 
algebra, one obtains what was claimed. 
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5.3. First Correction to the Variational Approximation 



The first correction to the variational procedure is given by the second cumulant ^\ The 
computations are essentially a copy-paste of those done in the eikonal representation. Espe- 

(5.42) 



cially, formula (4.78) still holds, 



A 



(3-1) 



\Xray) (Xray) 



The evaluation of (Xray) is done in appendix 



A. 



one has, for n = 2, 



exp 



2m 



Pi ■ P2 Fl - t: 



t2\±] ) - exp 



2m 



[pLI^iI + pLI^2|] 



With the definition (5.7|) of the effective potential it reads 

{Xray} 



j dti J dt2 J Ylj^K(u)iPi) ^^5^P' ■ 



exp 



2m 



Pi ■ P2[|^l - til - - \t2\±] 



(5.43) 



The only differences with respect to the eikonal representation are thus the presence of the 
effective potential and the terms \ti\± in the exponent. One sees that the first two terms in 



the developpment of this exponential give rise to X^^ squared and 2iX[^ respectively, so 
that the second cumulant is that expression, starting from the third term in the expansion of 
the exponential function. 



'(3-3) 



5.4. High Energy Expansion 

As in the eikonal representation, we fix the velocity k/m and expand the trajectory and the 
corresponding phases in inverse powers of k. Since we will not investigate in much detail 
beyond the first order, we will not make any difference between k and K. 
We first have a look to the properties of the effective potential in this approach. From its very 
definition (5.7), one has 

Kw(p) = Vip)exp (-z^|t|pl) ^ Vip) - z^|t|ply(p), (5.44) 

so that for large K, or for small times, the effective potential is V with a complex correction, 

V.it)-V + t^\t\A^V. (5.45) 



Especially, the imaginary part of the variational trajectory behaves, under such conditions (for 
example near the scattering region), under the influence not of V, but of the (perpendicular) 
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Laplacian of V. This imaginary trajectory is clearly a second order effect. Indeed, the zeroth 
order zo(t) is the usual straight line 



K 



q 



Zo{t) =h+-t + 

m 2m 



while the first order term zi(it) is readily obtained through 



zi W = j ds VV{zo{t)) [\t - s| - \t\ 



(5.46) 



(5.47) 



where the zeroth order term V of V^(s) had to be used. Thus zi(t) is still real, while in the 

second order term of z(t) would appear the first order term of V^(s), which is complex. 

Let us have a look at the scattering phases to first order. The factor 1/m = v/K in front of 

^(3-1) 

means that, up to higher order terms. 



(3-1) 



ds / dt VV{zo{s))VV{zo{t)) [\t - s\ - \t\± - |s| 



(5.48) 



(3—1) 

On the other hand, the phase Xq can be written as 



X, 



(3-1) 



dt 



dp 

{2ny 



V{p) exp -— pi|t| exp [ip ■ z{t)] 



2m 



(5.49) 



dt 



dp 



2m 



Vip) 1 - — pl|t| exp [^p ■ zo(t)] (1 - ^p ■ zi(t)) , 



(3—1) 

and contains two real parts, one of zeroth order, that we will call Xo \ and one of first order 
which is recognized, using (5.47), to be 

dt 



dp 

(2^ 



ip V{p) zi{t) exp {ip ■ zo(t)) 



V 

2K 



j ds j dt Vr(zo(s))V\/(zo(t)) [\t -s\- \t\ 



(5.50) 



However, Xq^ contains to first order also an imaginary part ^' , given by 



(3-1) 



(3-1) 



V 

2K 



dt \t\ 



dp 



3?! V"(p)exp(ip ■zo(t)) 



V 



(5.51) 



where Af, is the two dimensional Laplacian operator acting on b in zo(t). To sum up, one has 

(5.52) 



after adding (5.48), (5.50), (5.51), up to higher order terms. 



(3-1) 



with 



(3-1) 

Xo 



- / dtV{zo{t)), 



(5.53) 
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(3-1) 



and 



(3-1) 

Xi 



V 

AK 



ds I dt Vy(zo(s)) ■ Vy(zo(t)) [\t - s\- \t\^ - |s| 



(5.54) 



(5.55) 



These three phases are identical with those obtained in the high-energy hmit in [7]. We refer 
to it for a more complete discussion. 
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6. Numerical Implementation and Results 



In this section we test numerically the approximations of the T-Matrix we obtained in this 
work. More concretely, we compute numerically the differential cross section 



a{n) = \fm 



m 



IT- 



(6.1) 



where Ti^f are the expressions (3.18) and (3.19). 

These are tested against the eikonal expansion of Arbabanel and Itzykson, and partial wave 
calculations. A sufficient number of partial waves were considered, so that we could consider 
them as the exact results. 



6.1. Eikonal Representation 

This is done for the case, first, of a Gaussian potential of the form 

Vg{x) = Voexp {-ax"^) , 
and second, for a Woods-Saxon potential 



Vws{x) = Vq 



The different parameters were chosen as 



1 + exp 



x-R 



k 
a 



4 fm" 



10fm-\ Vn 



-0.2 fm" 



1 fm~\ i? = 2 fm and a = 0.2 fm. 



(6.2) 



(6.3) 



(6.4) 
(6.5) 



Using the relation He ~ 200 MeV fm~^, these correspond to a particle of mass 2 GeV, with an 
initial kinetic energy of 160 MeV, interacting through a potential well of depth -40 MeV. 
We are in a high energy situation in both cases, and these parameters are precisely those for 
which the eikonal expansion was judged in [4] as being unsatisfactory. 

For our calculation, the variational trajectory x(i(:) needs to be accurately found, for each 
scattering angle and b. We determined it through iteration. Starting from the eikonal reference 
trajectory 



Xo(t) 



K 



m 



a sufficient number of iterations 



x„(t) := 



b H 1 

m 



1 
2m 



ds VV(x„_i(s))|t - s| 



(6.6) 



(6.7) 



was performed. Since the potentials are spherically symmetric, the motion takes place in both 
cases in the plane determined by b and K. The x-axis was taken parallel to b and the z-axis 
parallel to K. Figure [2] shows the successive trajectories x„(t) in the case of the Gaussian 
potential, for the given parameters and a particular choice of 6 = 0.1 fm. Note how the z-scale 
was stretched to see the effect on the asymptotic momenta. Figure |4] shows the corresponding 
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trajectories for the Woods-Saxon potential. To evaluate numerically the integrals needed in 
this iteration procedure, we used a Gauss-Legendre rule with 32 points and a sufficient number 
of subdivisions of the integration interval. The trajectories were numerically stored at points 
ti, determined through the mapping 



TT TT 

2' 2 



where the '?/'j/(7r/2) are the Gauss-Legendre points, and 

m m 



tn 



k\/a 



or t 







k 



(6.9) 



is the typical time scale of the scattering process, for the Gaussian and Woods-Saxon potential 
respectively. 




Figure 2: The successive approximations x„(t) of the variational trajectory x(t), for a particu- 
lar value of 6 = 0.1 fm, and K = 2 fm~^. The potential is Gaussian, with the parameters given 
in the text. The z-scale has been largely stretched, so that the differences in the asymptotic 
momenta can be clearly seen. 

Gaussian Potential. With the help of the so-obtained trajectories, the phases Xq and Xi 



were evaluated according to equations (4.42) and (4.51), and the two-dimensional integration 
over the impact parameter b was reduced to a one-dimensional integration, using the fact that 
the phases depend only on b. All these integrations were performed using the Gauss-Legendre 
rule we mentioned above. 

The specific form of the Gaussian potential permits also the analytical computation of the 
first correction to the variational approximation without much trouble. We have already seen 
that the second cumulant takes the form 

\2 = {xk)-X'o-2zX,, (6.10) 
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while in appendix [C]] we evaluate, for the Gaussian potential, 



^11 

^11 



exp I -2z«7,|^^M^^M ) , (6.11) 



with 



11l(^-fel) = -|*.-fc|. (6.12) 

m 

However, the oscillatory character of this second cumulant proved to be very troublesome in 
trying to evaluate it numerically. With the Gauss-Legendre rule we used for the evaluation of 
the first cumulant, it turned out not to be possible to achieve in reasonable time its computa- 
tion to a sufficient accuracy. What we did to circumvent the problem was to use the numerical 
package DCUHRE [13], that uses an adaptive scheme for multi-dimensional integration, and 
permitted an accurate evaluation of the second cumulant for most values of b. For the largest 
values of b, we used a third method, that relies on the fact that the variational trajectory is 
essentially the same as the eikonal straight-line reference trajectory. More details can be found 



in appendix E 



Results are shown in figure |3| Presented in the plot is the exact cross section, evaluated 
with a sufficient number of partial waves, the eikonal (AI) approximation with the first correc- 
tion included, and the variational approximation, with and without the first correction. The 
upper panel shows, on a linear scale, the relative deviation with respect to the exact scattering 
amplitude as a function of angle. Although the variational approximation, without the first 
correction, suffers the same deviation as the eikonal approximation at about fifty degrees, it 
clearly does better at larger angles. Especially, when the first correction is included, the im- 
provement in this region over the eikonal approximation is impressive, as shown clearly by the 
upper panel. This is hkely to be traced back to the appearance of the real term in A2. 



Woods-Saxon Potential. We performed the same computations in the case of the Woods- 
Saxon potential, although we did not consider the first correction to the variational approxi- 
mation. We used only Gauss-Legendre integration for each needed integral, with a sufficient 
number of points and subdivision of the integration intervals. Results are presented in figure 

El 

Due to the edge in the Woods-Saxon potential, the cross section shows a diffraction structure, 
different from the Gaussian potential. While the eikonal approximation levels off and does 
not manage to reproduce the diffraction peaks at large angles, the variational approximation 
does it, although it gets more and more out of phase. The action of the first correction to the 
variational approximation would be interesting to study, albeit numerically quite involved. 
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20 40 60 80 100 120 

e (degrees) 



Figure 3: The differential cross section a as function of the scattering angle, for the Gaussian 
potential. The red line shows the exact, partial-wave results, the blue line the eikonal (AI) 
phase with the first correction. The other two lines show the variational approximation in the 
eikonal representation, without (magenta line) and with the first correction (green line). A 
sensible deviation is found at about 50 degrees, above which the variational approximation is 
then much closer to the exact result than the eikonal approximation. The upper panel shows 
the relative error in the scattering amplitude, with respect to the partial-wave calculations, 
for the eikonal approximation with first correction, and variational approximation with first 
correction (green line), for scattering angles up to sixty degrees. Note that at this value, the 
eikonal approximation is off by more than 10 percent, while the variational approximation still 
is within 2 percent. 
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Figure 4: The same as figure [2| in the case of the Woods-Saxon potential. The significance of 
the parameters i? = 2 fm (width of the potential well) and a = 0.2 fm (width of the edge of 
the well) can be clearly seen. 



6.2. Ray Representation 

We computed also the cross section in the ray representation in the case of the Gaussian 
potential. The variational trajectories were again obtained through iteration. Starting from 



h + —t 

m 



q 

2m 



t\ 



the successive approximations z„(t) were found by setting 

Z„(t) ■=h+ -t+^\t\- — ds VK(s)(Zn-l(s)) [\t -S\- \t\ - \sW 

m m 2m J 



(6.13) 



(6.14) 



The same Gauss-Legendre rule as in the eikonal representation was used to find these trajecto- 
ries. However, there are three significant differences with respect to the eikonal representation. 
First the motion is not in a plane anymore. Second, the trajectories are complex numbers, 
and third, one must perform numerically the full 2-dimensional integration over b. 
Also, the effective potential Vcr(t) is a new feature. For the Gaussian potential Vq, it is readily 
computed from equation (5.7): 



Vo 



dp 

(27r)3 
Vo 



7r\3/2 



a J 



) exp 



1 

4a 



p2 - J(T(t)pl + ■ z(t) 



2aa{t) 



exp 



-a 



2aa{t)_ 



(6.15) 



Different from the case of the eikonal representation, we performed all the other integrations 



with the help of the DCUHRE package. The first correction is also obtained in appendix E. 



Results are plotted in figure pi The variational approximation in the ray representation without 
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Figure 5: The differential cross section a as function of scattering angle in the case of the 
Woods-Saxon potential. Shown are the exact result given by the partial-wave expansion (red 
line), the eikonal (AI) phase with the first correction included, and the variational approxima- 
tion in the eikonal representation (first cumulant only). Again, the variational approximation 
shows improvement at large angles. While the eikonal approximation is clearly inconsistent 
for such angles, the variational approximation is able to keep track of the diffractive structure, 
although it gets slowly out of phase. The upper panel shows the relative error in the scattering 
amplitude / with respect to the partial-wave computations. 
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first correction (magenta line) shows, plausibly due to the presence of the real term in the 
scattering phases, a very similar behavior as its counterpart in the eikonal representation, 
but with the first correction included. We see our prediction that the ray representation will 
be more accurate clearly confirmed. With the first correction included, the result of the ray 
representation is even more impressive. As the upper panel shows, the approximation is within 
two percent of the exact result, all the way through ninety degrees, and only then begins to 
deviate. 




20 40 60 80 100 

e (degrees) 



Figure 6: The differential cross section a as function of scattering angle in the case of the 
Gaussian potential. Shown are the exact result given by the partial- wave expansion (red line), 
the eikonal (AI) phase with the first correction included, and the variational approximation 
in the ray representation, without (magenta line) and with first correction (green line). The 
upper panel shows the relative deviation in the scattering amplitude. With the first correction 
included, the variational approximation in the ray representation is rather accurate, even at 
large angles. 
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7. Summary and Outlook 

In this master thesis, we have investigated a new variational approach to potential scatter- 
ing, starting from two path-integral representations of the T-matrix, where one integrates 
functionally over velocities. Although it shares some common properties with the eikonal ap- 
proximation, it possesses new and interesting features, and has a greater domain of validity. 
The Feynman- Jensen variational principle we used, together with our choice of trial action 
(the free action plus a linear term to be determined), select the best reference trajectory the 
particle travels along, given in the form of an integral equation of motion. The T-matrix is 
then given by summing up phases along this path, as in the eikonal approximation. Unlike its 
eikonal counterpart, this trajectory is in the first representation a real scattering trajectory 
that obeys Newton's equation of motion. In the second representation, interestingly, the path 
is complex, and obeys a similar classical equation of motion, albeit under the influence of a 
new, efi^ective potential which has an explicit dependence on time. The high-energy behavior of 
the approximations was then investigated. Although the recovery of the leading and next-to- 
leading terms of the eikonal expansion was expected, the fact that the second order imaginary 
term was also found came as a pleasant surprise. Numerical tests of this approximation in 
a high energy situation have then shown that it does better at large angles than the eikonal 
approximation of Abarbanel and Itzykson. The numerical evaluation of the first correction to 
the variational result (the second cumulant) turned out to be very challenging. However, it 
became feasible with the use of adaptive integration procedures. This was really worth a try, 
since, especially in the second representation, the agreement between the exact results and our 
approach is remarkable. 

Naturally, this method also has its shortcomings. For example, the choice of trial action 
is restricted to those at most quadratic in the velocities. We studied in this work one simple 
choice of linear Ansatz. with the physical motivation that such a linear term in the T-matrix 
would emulate the structure of the high energy expansion, albeit without being subjected 
to high energy conditions. It would be interesting to study another, more general kind of 
quadratic trial action, where one also determines the quadratic part of the trial action from 
the variational principle. Inevitably, this will lead to more involved, coupled equations, that 
are more difficult to solve. Also, we have been unable to find a way to analyse the behavior of 
the approximation for low energies, although one may expect that the variational approxima- 
tion could also be meaningful in this kinematic regime. 

As further investigations of this formalism one could propose the extension to relativistic 
scattering, as well as applications in many-body processes. 
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Appendix 



A. Computation of (x 



^ray 



In this section we calculate the expectation value (Xray) of the phase xk in the ray represen- 
tation, that was used without proof in the text: 



i=l 



As in the text, whenever a _L symbol is found, like here as a companion of \ti\ and \tj\, it is 
meant that only the part that is perpendicular to K of the 3-dimensional vector that multiplies 
it is to be taken into account (i.e. Pi|tj|_L = Pj±|tj|). Similarly for the y symbol. We recall that 
in the ray representation, xk is given by 

XK(b, w,w) = - j dtv{h+ + x,(t) - x,x(t) - a;^(0)^ . (A.2) 

The computation is very similar as in the case of a 3-dimensional antivelocity w. We first start 
with a Fourier transformation of the potential, and solve the functional v and w integrations. 
Let us consider first the v integration. 

e^f (v>v)+»(B,v) ^.1 J2Pijds v(s) [sgn(ti - s) - sgn^(s)] j 

J ds |B(s) + ^^Pi[sgn(t,-s)-sgn^(-s)]| j . (A.3) 



2m 



= exp 

Similarly, the w integration results in 



'^''P ' 2^ / "^M ^^'^ " 2 > ^ P.||Sgn(-3) \ I . (A.4) 



As already encountered, the and again cancel the normalisation ttlq. Let us consider 
more closely the terms in the exponent that multiply pj ■ Pj and are free of B's and C's: 

j ds (sgn(ti - s)sgn{t,j - s) + U - sgn_L(-s) [sgn(tj - s) - sgn(ti - s)] - 1||) 
= -2{\ti-t,\-\ti\^-\tj\A_). (A.5) 



Here we have used the relation (2.38) only. It is thus seen that all divergent quantities cancel. 



We further define the trajectory z(t) as 

z(i) = b + Pr^^t ^J-[ds B(s) (sgn(t - s) - sgn ± (s)) + C(s)sgni| (s). (A.6) 
m 2m J " 
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Getting the factors right, and plugging into {Xray)i reduces to 

/n /■ " ^ / 

n / n T^^^p^) exp z ^ p, ■ z(to 
i=i "'1=1'^'' \ i=i 



exp 



i=l 

which was the formula to be proven. 

B. Reduction Formulae 

B.l Reduction of Xi(^) 

We prove in this section the formula for the first order phase Xi{^)i for spherically symmetric 
potentials, which was mentioned in the text: 



If . d 



00 







Since we are dealing with a spherically symmetric potential, we will always write, in such 
expressions, V{ri) instead of V^^Jb"^ + Zf), for the sake of clarity. To this aim, we start from 
the general formula (4.97), and effectuate the substitutions Zi = ^t, and Z2 = ^s. We then 
obtain 



Xiib) = ^ J dZ, J dZ2 V\/(ri) • Vy(r2)|t - s\. (B.2) 

Since b and K are perpendicular to each other, and K points along the third axis, it is 
immediately recognized that dsVlh + ZK.) = dzV{r). We furthermore choose b to point 
along the first axis, so that diV{r) = dbV{r) also holds. We will make extensive use in the 
following calculations of the relation ZdbV{r) = bdzV{r), which is a special case of Xidj = Xjdi, 
that holds for any spherically symmetric function. 
We thus have 

Xi{b) = ^ j dZ, j dZ^ {dz,V{r^)dz,V{r2) + dbV{r{)dbV{r2)) \Z^ - Z<,\. (B.3) 

We are allowed to integrate over Z\ and Z2 by parts in the first term, while for the second 
term we first write {Z\ — Z2) sgn(Zi — Z^) for \Z\ — Z2I, and use Zd}y{r) = hdzV{r). We 
obtain 

Xiih) = ^ j dZ^ i^-2V\r) + 2bJ dZ^ ^z,V{r^)^bV{r2)sgn{Z, - Z^)^ (B.4) 
where we have used 

dz,dz,\Z, - Z2I = -dz,sgn{Z, - Z2) = -25{Z, - Z^). (B.5) 



J. Carron : Variational Methods for Path Integral Scattering 



53 



The second term can now also be integrated by parts, and, after recognizing that 

2V{r)dbV{r) = dbV\r), (B.6) 

one finally has 

Mb) = ^,JdZ (-2F^(r) - 2b^^V\r)^ . (B.7) 

If we use the symmetry of the integrand to integrate from to oo, one gets exactly what was 
to be proven: 

Xiib) = --^ f 1 + 41 r dZ V\r). (B.8) 



db 



B.2 Reduction of cu 



(3-3) 



In this section we derive the reduced expression of oo^ '^\b), for spherically symmetric poten- 
tials, that was mentioned in the text : 



8K^ 



2Xo^ + x'o {xo + bx'o) + I 



XoXo 



(3—3) 

This is given by the first term in the expansion of the second cumulant, evaluated along 

the straight line trajectory. 



1 



dZi / dZ2{Zi- Z2y^^^JV{h + ZiK)^,^JV{h + Z2K), (B.9) 



where a sum over the i and j indices is implicit. 

We will use the same conventions as in appendix B.l The integral consists of a few terms 
that need to be calculated. Let us consider them one by one. 



case i = j = 2: We will call this term Ai. Since 



(B.IO) 



it is the only non-zero term where one of the two indices takes the value two. It reads then 



Ai{h) = J dZi J dZ2 {Zi - Z2 
= '2 j dZiZl 



^2 V\r,) V\r2) 
ri r2 



V\r,) f V'{r2) 



r2 



(B.ll) 
(B.12) 



We used in the second line the symmetry properties of the integrand. The Z2 integral is, up 
to a factor of b and v, the derivative of Xo{b), and the first can be integrated by parts with the 
help of Z V' /r = dzV, so that we obtain 



(B.13) 
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case i — j — b: This term, that we call A2, reads 

M ^ j dZi j dZ2{Zi-Z2)%dbV{n)dbdbV{r2) 
= 2 y dZi ZldhdhVin) j dZ2 dhdhV{r2). 



(B.14) 
(B.15) 



The second integral is recognized to be the second derivative of Xo: while the first, using 
ZdbdbV = (1 + bdi,)dzV, is integrated by parts. We get 



Mb) = vxoj dZ [l + bd,]V{r) 
= -2^^^Xo(Xo + &Xo)- 



(B.16) 
(B.17) 



case (i, j) = {b, Z) and {Z, b): These two terms are again readily computed with integration 
by parts : 

A^(b) ^ 2 j dZ, j dZ2{Z,-Z2ydbdbV{n)dbdhV{r2) (B.18) 

= -4(^1 dZZdzdbV{r)^ (B.19) 
= (B.20) 
Since any other term than those three, obviously vanishes, we readily have proven our claim : 

ujt'\b) = Ar{b) + A2{b) + A,{b) 



2xif + Xo (xo + ^'Xo) + ^XoXo 



(B.21) 



B.3 Reduction of ^^(&) 



(3—3) 

This section is devoted to the derivation of the representation of X2 ■> the case of a 
spherically symmetric potential function, that was presented in the text: 



xr'(6)= - 



1 / 5, d 1,, (P \ f°° j^,,,, , 



8K^ 



Xo(Xo)' + 3KX0)' 



(B.22) 



We will use throughout this section the same notation convention as in the last section. The 
expression to be reduced is 



^ J dZ,,2,3 didjV(n)diV(r2)djV(r3)\Z, - Z2\\Z, - Z^\. (B.23) 



It is seen immediately that no terms with an index i or j — 2 contributes to the sum. Thus, 
there are three terms left to be computed. 
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case i — j — Z: This term, that we define as Ai{b), is the less involved. We integrate by 
parts to obtain : 

A,{b) := j dZ,,2,3dz,dz,V{n)dz,V{r2)dz,V{rs)\Z,-Z2\\Z,-Zs\ 

= J dZ,^2,3 dz,dz,V{n)V{r2)V{r2)sgn{Z, - Z2)sgn{Z, - Z^) 

= -2y" dZ^,2^z,Vin)Vir2)V{r^)sgniZ^-Z2) 

There we only made extensive use of the rule \x\' — sgn(x) and sgn(a;)' = 25{x). A last 
integration by parts permits us to conclude 

Ai{b) = A j dZV^r). (B.24) 

case = {b, Z) and (Z, b): These two terms consist of 

A2{b) := 2 J rfZi,2,3 dzAV{n)dz,V{r2)di>V{rs)\Zi - Z2\\Z^ - Z^\. (B.25) 

There we will encounter a function that plays an important role. We define it by the following 
equation, 

X{Z, b) := j dZi sgn(Z - Zi)V{ri). (B.26) 

It has two nice properties, which we will use many times: 

X{±oo, b) = ^vxoib) and dzX{Z, b) = 2V{r{Z, b)). (B.27) 

From now on we will not write explicitly the b-dependence of X, and write a dot for d^. 
After one integration by parts over Z2, one obtains 

A2{b) = 2J dZ,,sX{Z,)Vin)Vir3)dz,V{n)\Z^-Zs\ 

- -^j dZ X{Z)X{Z)V{r)-2 J dZi,sV{n)V{rs)V{ri)\Zi - Z^\. 

We leave the term on the left as it is for now. The second term can be further reduced. To 
this aim, we decompose \Zi — Z^l in sgn(Zi — Z^^iZx — Z3), and we effectuate the integration 
using Zdjj — bdz- After some algebra, one obtains finally 

A2{b) = J dZ V'{r) - 2F{b), (B.28) 

with 

F{b) :=- [ dZ X{Z)X{Z)V{r). (B.29) 
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case i = j = b: This quantity reads 

As{b):= J dZi,2,3\Zi-Zs\\Zi-Z2\V{ri)V{r2)V{rs). (B.30) 

First we consider the Zi and Z3 integrations. We have 

j dZ2V{r2)\Zr- Z2\ = j dZ2V{r2)sgn{Z^- Z2){Z^- Z2) (B.31) 

= ZiX(Zi) - 26F(ri), (B.32) 

so that ^3 can be rewritten as 

^3(6) = j dZV{r){zX{Z) -2hV{r)y (B.33) 

Further algebra involving integration by parts then reduces ^3 to 

A,{h) = J dZ V^r) + lb'^J dZ V^r) + 4G(6) + Hib), (B.34) 

where G and H are defined as 

G{b) = b'^ j dZV{r)X{r)dzV{r) (B.35) 

H{h) = j dZ Z{l + bdb){dzV{r))X\Z). (B.36) 

The last step is to recognize that 4G — 2F + if is none other than 

AG-2F + H = xoix'of + KXof + V\r). (B.37) 

This is done of course also by integrating by parts, but this time the integrand evaluated at 
infinity contributes, making ^iid its derivative to appear. By summing over the different 
quantities so computed, one obtains exactly what was claimed. 

C. The Second Cumulant for a Gaussian Potential. 

In this section the computation of (x^) is done in the case of a Gaussian potential of the form 

y(x) = ^oexp (-ax^) . (C.l) 
The second cumulant is then, in both the eikonal and ray representations, given by subtracting 



the phases Xq and 2iXi in their respective representations (see equations (4.78) and (5.42)). 



After plugging this form of the potential into the formulas for (x^) we already derived, we get 

(^K> = < i'-yj j'^'^ji^j^. xK(i„i.,p„p.) (C.2) 
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in both cases, with 



and 



= exp (P? + P2) + ^Pi ■ + ip2 ■ z{t2) + ^Pi ■ P2|ti - i2|| , 

XK"^^(^l,i2,Pl,P2) = XK'^\tl,t2,Pl,P2) ■ 

■ exp (pLI^iI + pLI^2|) - ^Pi± ■ P2± (|ti| + |t2|)| . 

In these expressions, z(t) means the reference trajectory that we denoted x(t) in the eikonal 
representation, or by the same letter z in the (3-1) case. These are 6-dimensional Gaussian 
integrals and can be thus given in closed form. In order to effectuate this more efficiently 
in both cases, we first consider the integral over one component of pi and the corresponding 
component of p2- These are all integrals of the form 

/ ^ y ^ exp {-aipj - a2pl + ipiZi + ip2Z2 + hpiP2) , (C.3) 

with different coefficients a and h for the perpendicular and parallel components. These 2- 
dimensional integrals are easily done and are equal to 



1 / a2zl + aizl + hziZ2\ I 1 ,^ 
2;^^"H 4a,a.-6^ J V i^^^^^' ^^''^ 

Eikonal Representation: The eikonal representation does not distinguish between parallel 
and perpendicular components, so that for each pair of integrals, 

1 i 

ai = 02 = — — and h= — \ti — t2\. (C.5) 
4a 2m 

Besides, we define a new quantity 7|j as being 

a 

7||(ti,t2) := -|ti-t2 , (C.6) 
m 



such that 4aia2 — 6^ is then simply (1 +7||)/4a^. Equation (C.4), applied to each component 
of pi and p2, then says 

3/2 / 

7?{ti)+7?{t2) 



(Xk) = j dti j dt2 { ) ( ~" 



l + 7||y V l + 7j| 

■exp|-2.a7„^^f^'l. (C.7) 
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Ray Representation: In this case, different coefficients a and b correspond to the parallel 
and perpendicular components. The integration over the parallel components is the same as 
in the eikonal representation, but for the perpendicular components, we have 



1 



4a 2m 
With the introduction of 



It I a - ^ + ^ 
4a 2m 



l'^2|, b = — {\ti - t2\ - \ti\ - \t2\) . 

2m 



7i,2(i 



a 



m 



we have 



4aia2 — b'^ 



4a^ 



and 73(^1, i^2) := 7|1 - 7i - 72, 



[(l + 2z7i) (1 + 2^72) +73] 



and finally, according to equation (|C.4|) 



' ray 



dt') 



l + 7f 



[(l + 227i)(l + 2z72)+7|] 



21-1 



• exp —a 



exp ( — 2za7; 



zi(ti)(l + 2z72) + zi(t2)(l + 2i7i 



(l + 227i)(l + 2z72)+7| 

Z±(ti) ■ Z_L(t2) 

7| 



exp —a 



1 + 2z7i) (1 + 2^72 

Z2(ti)+Z2(t2 



l+7f 



exp — 2za7|| 



Z||(ti) • Z||(t2) 

l + 7f 



(C.8) 

(C.9) 
(C.IO) 



(C.ll) 



D. Angle Corrections to the 2nd Order Scattering Phases 

The aim of this section is to derive two results that were mentioned in the text that contributed 
to the second order scattering phases U2 and T2 in a high-fc expansion of the scattering matrix 
in the eikonal representation. The phases X2^ (^) ^md 002^ (b) were defined as the imaginary and 
real part of the second order term arising because of the presence of the cos(^/2) term in the 



expression (2.32) for the T-matrix in the eikonal representation. Similarly, the imaginary and 
real part of the second order term arising because of the presence of the very same term in 
the leading phase were named X2^{b) and (^2^(6). The identities proven in this section are 



and 



^X2^ - ^2' = ^ i^^bXo - (Vao)') (D.l) 
ZX2' - 00^2' = ^ (2 (Vao)' + XoAao + ^((Vao)\o - A^o)) , (D.2) 



8F 

where Xo denotes the standard eikonal phase 



m 

Xoib) = -j I dZV{r{b,Z)). (D.3) 
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We first prove (D.l ). To this aim, we first recognize that, up to second order in 1/k, cos(^/2) = 
(1 — g^/S/c^). This permits us to rewrite the T-matrix as (in the subsequent development an 
equal sign always means equahty up to terms of order four and higher.) 

,2 



7^(3-3) _ • 



k 



m 



d'b 1- 



8P 



(D.4) 



Since we are interested in the second order terms, it is sufficient to consider only the leading 
phase Xo the exponential function. The factor can be trivially converted to a Laplacian 
operator, and then brought on the other side by an integration by parts. 



7^(3-3) _ • 



k 



m 



Since the identity 



holds, one obtains 



T, 



(3-3) 



m 
. k 

^ 

m 



1 

8F 



(D.5) 



(D.6) 



(D.7) 



J d% e^'^ '' exp jzxo + ^ [^A^o - (V^o)'] } - 1 



which is the first of the expected results. To prove (D.2), the procedure is very similar. From 

(D.8) 



(D.3) it follows that 



which in turn implies 



(3-3) 



K 

m 
K 
m 



Xo = Xo 1 1 + ^ 



J d'be-^^ '' (exp|zxo( 



1 + 



1 



d^b e-'"'-^ I e^^« I 1 + 



IP 



1 



(D.9) 
(D.IO) 



Again, equalities are to be understood to hold up to fourth order. We now use the same trick 
as we just did, i.e. we convert the g^-term to a Laplacian operator and integrate by parts: 



7^(3-3) _ • 



K 



m 



d% e-^'i-^ <^ e^^° I 1 



A small calculation shows 

Afe [xo e*^°] = e^^° (- (Vac)' Xo + A^o + i (2 (V^Xo)' + XoAao)) , 
and subsequently. 



(D.ll) 



(D.12) 



7-(3-3) = i— f d'^be'''^'^ ^g(iXo-^[-(V6Xo)\o+A6Xo+i(2(V6Xo)^+XoAbXo)]) _ ^ (D.13) 



which is what was to be proven. 
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E. Numerical Evaluation of the Second Cumulant 

In this section, we briefly explain the difficulties encountered in evaluating numerically the 
second cumulant for the variational approximation in the ray representation for the Gaussian 



potential (6.2). As we already said in the section on numerical results, we were unable to eval- 
uate the main component {xk) of the second cumulant accurately, through Gauss-Legendre 
integration. Since the Gauss-Legendre points enter our program in every integration that had 
to be computed, the requested massive increase of the number of points would be too costly 
in time to be performed. Instead, we chose to adopt another integration scheme, the adap- 
tive integration package DCUHRE, to calculate the second cumulant. Since the variational 
trajectories were stored only at places corresponding to the Gauss-Legendre points, a linear 
interpolation between these points was also done. As a typical example of the problem en- 
countered, figure [7] shows the modulus of the real part of (xic) computed in these two different 
ways (blue and red curves), for a special value oi K = 2 fm~^, and for integration parameters 
adapted to our situation, as a function of b. It shows that two orders of magnitude in precision 
were gained through this procedure. 

However, this was still not enough, since the adaptive integration also levels off after a given 
value of b, that was still to be accounted for in the computation of the differential cross section. 
To solve this problem we used the fact the trajectory is sufficiently close to a straight line to 



approximate (xk) by its value given by (6.11), with 



h + —t instead of x(t). (E.l) 

m 

In that case, after passing to the variables 

T:=^{ti+t2), and t ■= ti - t2, (E.2) 

the integral over T is a pure Gaussian and can be analytically solved. The 1-dimensional 
integration that is left was then numerically evaluated and compared to the full 2-dimensional 
integration involving the variational trajectory. It leads to the magenta line in figure [7| As 
one can see, the agreement between this curve and the one representing the two-dimensional 
adaptive integration is good enough, before the adaptive integration fails. 
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Gauss-Legendre Integration 
Adaptive Integration 
1d Al 




b (fm) 



Figure 7: The modulus of the real part of (xic) for the Gaussian potential, as a function of b, 
evaluated numerically by Gauss-Legendre integration and adaptive integration, with optimal 
control and accuracy parameters. The third, magenta line, is the 1-dimensional version, as 
explained in the text. Both the Gauss-Legendre as well as the adaptive integration proved 
unable to perform the integration up to sufficiently large values of b, although the adaptive 
integration levels off two orders of magnitudes later. 
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